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Fast  Approximation  Schemes 
For  Multi- Criteria  Flow, 

Knapsack,  and  Scheduling  Problems 

Hershel  M.  Safer  and  James  B.  Orlin 

Abstract 

The  solution  to  an  instance  of  the  standard  Shortest  Path  problem  is  a  single  shortest  route  in  a  directed 
graph.  Suppose,  however,  that  each  arc  has  both  a  distance  and  a  cost,  and  that  one  would  like  to  find  a 
route  that  is  both  short  and  inexpensive.  A  solution  in  this  case  typically  consists  of  multiple  routes  because, 
in  general,  no  point  is  both  shortest  and  cheapest.  Finding  the  efficient  frontier  or  tradeoff  curve  for  such  a 
multi-criteria  problem  is  often  difficult,  in  part  because  the  solution  set  can  have  exponentially  many  points. 
As  a  result,  multi-criteria  problems  are  often  solved  by  approximating  the  efficient  frontier. 

A  fast  approximation  scheme  (FAS)  for  a  multi-criteria  problem  is  an  algorithm  that  can  quickly  find  an 
arbitrarily-good  approximation  to  the  problem’s  efficient  frontier.  Necessary  and  sufficient  conditions  for 
the  existence  of  an  FAS  for  a  problem  were  introduced  in  [S094].  The  conditions  are  stated  in  terms  of  the 
existence  of  a  V-pseudo-polynomial  (VPP)  time  algorithm  for  the  problem.  A  new  form  of  reducibility  is 
also  introduced  and  shown  to  be  useful  for  studying  the  existence  of  FAS’s. 

This  paper  extends  the  work  of  [S094]  by  applying  it  to  a  variety  of  discrete  optimization  problems.  The 
Source-to-Tree  (STT)  Network  Flow  problem  is  introduced.  This  problem  is  interesting  because  it  generalizes 
many  commonly-treated  combinatorial  optimization  problems.  A  VPP  time  algorithm  is  demonstrated  for 
a  particular  STT  problem  and  is  used  to  show  that  the  problem  has  an  FAS. 

Results  are  also  derived  about  the  computational  complexity  of  finding  approximate  solutions  to  several 
multi- criteria  knapsack,  scheduling,  and  production  planning  problems.  Many  theorems  extend  previously- 
known  results  to  multi-criteria  problems.  For  some  problems,  results  about  problems  with  binary  variables 
are  extended  to  problems  with  general  integer  variables.  These  and  other  results  are  of  interest  even  for 
problems  with  only  a  single  criterion. 
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Glossary 

This  section  contains  brief  explanations  of  many  symbols  and  terms  that  are  used  in  this  paper  and  its 
companion  paper,  along  with  references  to  their  full  definitions. 

Symbols 

AV  (w),  ATI  (w),  AM  (w)  The  set  of  additively  separable  functions  that  are,  respectively,  in  P  (w),  7?.  (w), 
and  M{w).  ([2],  section  3.2) 

C  {Al/CG/MG/NM/FS)  ([2],  section  2) 

C  The  subclass  of  C  consisting  of  those  networks  in  canonical  form.  ([2],  section  4.1) 

Dui  (u,  v)  u  dominates  v  w.r.t.  w.  ([1],  section  3.3) 

(u,  v)  u  e-dominates  v  w.r.t.  w.  ([1],  section  3.4) 
ej;  A  vector  of  length  k  in  which  each  component  is  one.  ([1],  section  3) 

(a;)  The  excess  flow  into  node  v  that  arises  from  the  flow  x.  ([2],  section  2) 

£  A  non-negative  bound  on  relative  error.  ([1],  section  3.4) 

=  vpp  VPP  equivalent.  ([1],  section  5.1) 

fi  The  component  of  the  multi-criteria  objective  function  /.  ([1],  section  3.1) 

T  A  family  of  r-criteria  objective  functions.  ([1],  section  3.1) 

[//vj  The  function  /  scaled  by  v.  ([1],  section  4.2) 

G)-^)  The  directed  graph  defined  by  the  partial  order  ([2],  section  8) 

I  An  instance  of  a  problem.  ([1],  section  3.1) 

(5,  /,  u)  An  optimization  instance. 

{S,  f,u>,M)  A  feasibility  instance. 

L  (I)  The  length  of  instance  /.  ([1],  section  3.2) 
log(a;)  Base  2  logarithm  of  x.  ([1],  section  3) 

min{/,  M}  The  function  /  with  box  constraints  M.  ([1],  section  4.2) 

My  (I)  The  largest  value  of  instance  I.  ([1],  section  3.2) 

My  (I)  A  close  upper  bound  on  My  (7)  that  can  be  found  quickly.  ([1],  section  3.2) 
n  The  number  of  variables  in  a  problem  instance.  ([1],  section  3.1) 

I I  a;  1 1  The  infinity  norm  of  the  vector  x.  ([1],  section  3) 

w  A  direction  vector.  ([1],  section  3.1) 

u!^  The  dual  of  the  direction  vector  w,  obtained  by  reversing  the  direction  of  each  component  of  w. 

([2],  section  3.3) 

fl  The  set  of  ail  direction  vectors.  ([1],  section  3.1) 

fl'’  The  set  of  all  direction  vectors  of  dimension  r.  ([1],  section  3.1) 

'P  (ui),  TZ  (w),  M  (w)  The  set  of  functions  that  are,  respectively,  order-preserving,  order-reversing,  and  order- 
monotone  w.r.t.  u.  ([2],  section  3.1) 

n  A  problem.  ([1],  section  3.1) 

(^,  y,  w.  Opt)  An  optimization  problem. 

{'9,  F,ui,  Feas)  A  feasibility  problem. 

11/3 The  restriction  of  problem  11  to  instances  in  which  the  integers  in  are  bounded  by  the  polynomial 
function  p(  )  of  the  instance  length.  ([1],  section  5.3) 
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r  The  number  of  criteria  in  a  problem.  ([1],  section  3.1) 

ocvpp  VPP  reduces.  ([1],  section  5.1) 

A  family  of  feasible  sets.  ([1],  section  3.1) 

S  The  feasible  set  of  a  problem  instance.  ([1],  section  3.1) 

gcov  ^  i^cov  covering  forms,  respectively,  of  the  feasible  set  S  and  the  family  $  of  feasible  sets.  ([2], 
section  6) 

'  !  '  !  '  /■)  Specification  of  a  class  of  STT  networks.  ([2],  section  1) 
u{S)  A  vector  of  upper  bounds  on  the  components  of  points  in  the  feasible  set  S.  ([1],  section  3.1) 
Umaa:(S')  The  largest  component  of  u{S).  ([1],  section  3.1) 

.2+  The  set  of  non-negative  integers.  ([1],  section  3.1) 

In  an  STT  network  G  =  (A,  A,  c,  m,  6):  ([2],  section  2) 

N  The  set  of  nodes. 

A  The  set  of  arcs, 

c  The  vector  of  arc  capacities, 

m  The  vector  of  arc  multipliers. 

6  The  vector  of  node  demands. 

In  the  construction  of  the  first  VPP  algorithm;  ([2],  section  4.1) 

T  The  tree  corresponding  to  the  graph  G. 

p  The  root  of  the  tree  T. 

I  (w)  The  postorder  label  of  node  v. 
dy  The  number  of  children  of  node  v  in  T. 

T[v,j]  The  subtree  of  T  defined  by  v,  the  first  j  children  of  v,  and  all  the  descendants  of  those 

children. 

G[v,j]  The  subgraph  of  G  defined  by  node  0  and  the  nodes  of  G  that  correspond  to  the  nodes 
of  r[w,  j]. 

j,  k)  The  maximum  excess  flow  into  node  v  over  all  flows  x  that  are  feasible  w.r.t.  G[v,  j]  and 
have  /  {x)  —  k. 

Terms 

Additively  separable  Describes  a  function  f  :  S  2’’+,  where  S  C  Z"+,  that  can  be  written  as  f  (x)  = 
E"=i  fj  (®l)-  ([2],  section  3.2) 

Arborescent  Describes  a  collection  of  sets,  any  two  of  which  either  are  disjoint  or  for  which  one  is  properly 
contained  in  the  other.  ([2],  section  7.1) 

/3-strongly  NP-hard  Describes  a  problem  11  for  which,  for  some  polynomial  p(  ),  the  problem  Up^p  is 
vVP-hard.  ([1],  section  5.3) 

Binary  family  A  family  of  feasible  sets  with  domain  {0, 1}.  ([1],  section  3.1) 

Box  constraint  Upper  bound  imposed  on  the  value  of  a  function.  ([1],  section  4.2) 

Box  constraint  neighbor  A  function  that  approximates  the  box-constrained  value  of  another  function. 
([1],  section  4.2) 

Closed  under  box  constraints  Describes  a  family  of  r-criteria  functions  for  which  applying  box  con¬ 
straints  to  any  function  in  the  family  yields  another  function  in  the  family.  ([1],  section  4.2) 

Closed  under  scrding  Describes  a  family  of  r-criteria  functions  for  which  scaling  any  function  in  the  family 
yields  another  function  in  the  family.  ([1],  section  4.2) 
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Domain  With  reference  to  a  family  ^  of  feasible  sets,  the  set  of  integers  from  zero  through  largest  value  of 
'U'maxiS)  over  sets  S  ([1],  section  3.1) 

Dominate  A  vector  u  dominates  v  w.r.t.  w  if  u  is  component-by-component  better  than  v  in  the  direction 
specified  by  w.  Written  (u,v).  ([1],  section  3.3) 

Efficient  Set  A  set  of  points  whose  values  define  the  efficient  frontier.  ([1],  section  3.3) 

£-dominate  A  vector  u  e-dominates  v  w.r.t.  w  if  each  component  of  u  is  no  more  than  a  factor  of  e  worse 
than  the  corresponding  component  of  v  in  the  direction  specified  by  w.  Written  Du^(u,v).  ([1], 
section  3.4) 

e-efficient  set  A  set  of  points  whose  values  are  within  a  factor  of  e  of  each  point  on  the  efficient  frontier. 
([1],  section  3.4) 

e-efficient  solution  A  minimal  e-efficient  set  for  an  instance.  ([1],  section  3.4) 

Exact  solution  For  an  optimization  instance,  a  minimum  cardinality  set  of  feasible  points  whose  values 
constitute  the  efficient  frontier;  for  a  feasibility  instance,  a  point  whose  value  achieves  the  target.  ([1], 
section  3.3) 

FAS  Fast  approximation  scheme.  Describes  an  algorithm  for  a  problem  11  which  for  any  e  >  0  and  any 
instance  I  €  H,  finds  an  e-efficient  solution  for  I  in  time  O  ({L{I)  /e)^^ ,  for  some  k  >  0.  ([1], 
section  3.4) 

Feasible  flow  A  flow  that  satisfies  the  flow  requirements  and  capacity  constraints.  ([2],  section  2) 

Feasible  set  A  bounded  set  S  C  .  ([1],  section  3.1) 

In-tree  partial  order  A  tree  partial  order  for  which  G(~>)  contains  a  node  toward  which  all  the  arcs 
of  G(^)  point.  ([2],  section  8) 

Largest  value  of  an  instance  For  an  optimization  instance,  the  largest  component  of  the  objective  func¬ 
tion  on  a  particular  box-constrained  superset  of  the  feasible  region;  for  a  feasibility  instance,  the  largest 
component  of  the  upper  bound.  Written  My  (/).  ([1],  section  3.2) 

Length  of  an  instance  The  number  of  bits  needed  to  represent  the  feasible  set  and  the  largest  value  of 
the  instance.  Written  L(I).  ([1],  section  3.2) 

Objective  function  A  function  f  :  S  Z’"'^ .  ([1],  section  3.1) 

Order-preserving,  order-reversing,  order-monotone  Describes  a  multi-criteria  objective  function  for 
which  increasing  some  components  of  its  argument  yields  changes  in  values  that  are  consistent  with  a 
specified  direction  vector.  ([2],  section  3.1) 

Out-tree  partial  order  A  tree  partial  order  for  which  G(-^)  contains  a  node  away  from  which  all  the 
arcs  of  G{'^)  point.  ([2],  section  8) 

Problem  A  collection  of  instances  of  the  same  kind  (i.e.,  optimization  or  feasibility)  with  identical  direction 
vectors,  feasible  sets  from  the  same  families,  and  objective  functions  from  the  same  families.  ([1], 
section  3.1) 

Quasi-closure  under  box  constraints  A  weaker  variant  of  closure  under  box  constraints.  ([1],  sec¬ 
tion  4.2) 

Quasi-closure  under  scaling  A  weaker  variant  of  closure  under  scaling.  ([1],  section  4.2) 

Scaling  Truncating  the  low-order  bits  of  a  function  value.  ([1],  section  4.2) 

Scaling  neighbor  A  function  that  approximates  the  scaled  value  of  another  function.  ([1],  section  4.2) 
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STT  Source-to-Tree.  ([2],  section  2) 

STT  network  A  generalized  network  which  has  a  particular  forest-like  structure.  ([2],  section  2) 

Sufficient,  deficient,  exact  flow  For  a  node,  refers  to  the  excess  flow  being,  respectively,  non-negative, 
non-positive,  or  zero;  for  a  network,  refers  to  the  flow  into  each  node  having  the  corresponding  property. 
([2],  section  2) 

Tree  partial  order  A  partial  order  for  which  the  undirected  version  of  G('^)  is  a  tree.  ([2],  section  8) 

VPP  V-pseudo-polynomial.  Describes  an  algorithm  for  a  problem  11  which  for  any  instance  /  G  11,  finds  a 
solution  for  I  in  time  O  (I)  My  ,  for  some  k  >  0.  ([1],  section  3.3) 

VPP  equivalent  Describes  two  problems,  each  of  which  is  VPP  reducible  to  the  other.  ([1],  section  5.1) 

VPP  reduction  A  VPP  algorithm  for  a  problem  that  uses,  as  a  subroutine,  a  VPP  algorithm  for  another 
problem.  ([1],  section  5.1) 

w.r.t.  With  respect  to.  ([1],  section  3.3) 

References  for  Glossary 

[1]  Hershel  M.  Safer  and  James  B.  Orlin,  Fast  Approximation  Schemes  For  Multi-Criteria  Combinatorial 
Optimization,  1994. 

[2]  Hershel  M.  Safer  and  James  B.  Orlin,  Fast  Approximation  Schemes  For  Multi-Criteria  Flow,  Knapsack, 
and  Scheduling  Problems,  1994. 
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1  Introduction 

The  solution  to  an  instance  of  the  standard  Shortest  Path  problem  is  a  single  shortest  route  in  a  directed 
graph.  If,  however,  each  arc  has  both  a  distance  and  a  cost,  then  one  might  ask  for  routes  that  are  simul¬ 
taneously  short  and  inexpensive.  A  solution  in  this  case  typically  consists  of  multiple  routes  because,  in 
general,  no  point  is  optimal  for  both  criteria.  Finding  the  efficient  frontier  for  such  a  multi-criteria  problem 
is  often  difficult,  in  part  because  it  can  have  exponentially  many  points.  As  a  result,  such  problems  are  often 
solved  by  approximating  the  efficient  frontier. 

Necessary  and  sufficient  conditions  for  the  existence  of  a  fast  approximation  scheme  (FAS)  for  such  a  problem 
were  introduced  in  [S094].  The  conditions  are  stated  in  terms  of  the  existence  of  a  V-pseudo-polynomial 
(VPP)  time  algorithm  for  the  problem.  In  addition,  a  new  form  of  reducibility  is  introduced  and  shown  to 
be  useful  for  studying  the  existence  of  FAS’s.  The  present  paper  applies  the  results  of  [S094]  to  analyze  the 
complexity  of  solving  a  variety  of  combinatorial  optimization  problems. 

The  Source-to-Tree  (STT)  Network  Flow  problem  is  introduced  in  Section  2.  To  the  knowledge  of  the  authors 
this  problem  has  not  been  considered  elsewhere.  This  problem  is  interesting,  however,  because  later  in  this 
paper  it  is  shown  to  generalize  a  wide  variety  of  commonly  treated  discrete  optimization  problems.  Several 
useful  classes  of  objective  functions  are  defined  in  Section  3,  and  a  VPP  algorithm  is  demonstrated  for  a 
class  of  STT  Network  Flow  problems  in  Section  4.  The  existence  of  this  algorithm  is  used  to  show  that  the 
problem  also  has  an  FAS.  The  boundary  between  easy  and  difficult  classes  of  STT  Network  Flow  problems 
is  explored  in  Section  5. 

The  STT  Network  Flow  problem  is  used  to  assess  the  complexity  of  finding  approximate  solutions  to  sev¬ 
eral  knapsack,  scheduling,  and  production  planning  problems  in  Sections  6  through  9.  Most  problems 
discussed  in  those  sections  have  been  previously  considered  with  only  a  single  criterion;  this  paper  extends 
the  previously-known  results  to  multi- criteria  versions  of  the  problems.  In  addition,  the  algorithms  described 
here  accommodate  the  general  integer  case  of  some  problems  for  which  only  the  0-1  case  has  been  considered 
until  now.  FAS’s  are  also  demonstrated  for  several  problems  that  have  not  been  considered  elsewhere. 

The  existence  of  an  FAS  for  a  general  lot-sizing  problem  is  demonstrated  in  Section  9.  Although  pseudo¬ 
polynomial  time  algorithms  were  previously  known  for  this  problem,  they  cannot  be  easily  transformed  into 
FAS’s.  The  results  of  Section  9  indicate  that  this  problem  has  a  pseudo-polynomial  time  algorithm  that  may 
be  more  useful  in  some  situations  than  are  previously-known  algorithms. 


2  The  Source- To-Tree  Network  Flow  Feasible  Set 

The  Source-to-Tree  Network  Flow  problem  is  introduced  in  this  section.  Not  only  is  this  problem  interesting 
for  its  own  sake,  but  it  is  also  useful  for  proving  the  existence  of  an  FAS  for  a  variety  of  other  problems.  As 
shown  in  Sections  6  through  9,  it  generalizes  quite  a  few  knapsack,  job  scheduling,  and  lot-sizing  problems. 
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Because  it  unifies  a  large  number  of  interesting  problems,  results  for  the  Source-to-Tree  Network  Flow 
problem  have  wide  applicability. 

This  section  explains  how  to  specify  feasible  sets  for  particular  classes  of  Source-To-Tree  Network  Flow 
problems.  Objective  functions  are  discussed  in  the  next  section.  A  VPP  algorithm  for  one  class  of  problems 
is  presented  in  Section  4.1;  results  for  other  classes  are  obtained  in  Section  5. 

To  start,  recall  two  definitions.  A  generalized  network  is  a  network  in  which  the  flow  that  enters  at  the  tail 
of  an  arc  is  multiplied  by  an  arc-specific  constant  before  it  emerges  at  the  head  of  the  arc[AM093].  Two 
directed  arcs  are  anti-parallel  if  they  connect  the  same  pair  of  nodes  but  are  oriented  in  opposite  directions. 

Intuitively,  a  source-to-tree  (STT)  network  is  an  m  -|- 1  node  generalized  network  that  satisfies  the  following 
conditions.  The  subgraph  induced  by  nodes  {1, . . . ,  m)  is  referred  to  as  the  tree  part  of  the  network.  The 
tree  part  of  the  network  has  a  special  forest-like  structure:  if  each  set  of  parallel  and  anti-parallel  arcs  is 
combined  into  a  single  undirected  edge,  the  result  is  a  forest.  Node  0  is  a  source/sink  node  and  is  adjacent 
to  some  other  nodes  in  such  a  way  that  the  entire  graph  is  connected.  Each  node,  except  the  source/sink 
node,  has  a  demand  or  supply. 

More  formally,  an  STT  network  is  defined  by  a  5-tuple  G  =  {N,A,u,p,b).  The  set  of  nodes  is  N  = 
{0,1,...,  m}  and  the  set  of  arcs  is  A  C  N  x  N.  Arcs  are  directed,  and  multiple  parallel  and  anti-parallel  arcs 
are  permitted.  If  only  a  single  arc  goes  from  node  v  to  node  w,  however,  it  may  be  represented  as  {v,  w)  for 
convenience.  The  network  is  connected,  and  the  arcs  in  A  have  a  particular  structure:  if  the  arc  directions 
are  removed  from  the  tree  part  of  the  network,  and  each  resulting  set  of  parallel  edges  is  replaced  by  a  single 
edge,  the  outcome  is  a  forest. 

The  vector  u  contains  the  flow  capacity  Ua  €  2'^  of  each  arc  a  E  A,  that  is,  the  upper  bound  on  the  flow 
allowed  on  arc  a.  The  vector  p  lists  the  corresponding  flow  multipliers  pa  G  2'^;  if  flow  x  enters  the  tail  of 
arc  a,  then  flow  paX  emerges  at  the  head  of  the  arc.  The  vector  b  specifies  the  demand  b^  E.  2  at  each  node 
V  E  N  —  {0}.  If  >  0,  then  «  is  a  demand  node;  if  <  0,  then  u  is  a  supply  node;  otherwise  b^  =  0,  and  v 
is  a  transshipment  node. 

Let  X  E  be  a  flow  in  G.  The  excess  How  into  node  v,  for  v  E  N  —  {0},  is  the  amount  by  which  the  flow 

into  V  exceeds  the  sum  of  the  demand  at  v  and  the  flow  out  of  v,  and  can  be  written  as 

Cy  (2:)  —  ^  ^  Pa^a  ^  by  . 

d^A'.  d^A". 

a=('w,v)  a={v,w) 

When  discussing  the  excess  flow  into  each  node  of  a  set  that  includes  node  0,  node  0  is  implicitly  excluded 
from  consideration  unless  otherwise  indicated.  The  flow  into  a  node  v  is  called 

sufficient  if  (x)  >0, 

deficient  if  Cy  (x)  <  0,  or 

exact  if  Cy  (x)  =  0  . 
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A  flow  X  is  said  to  be  suflBcient,  deficient,  or  exact  if  the  flow  into  each  node  v  £  N—{0}  has  the  corresponding 
attribute. 

Many  classes  of  STT  networks  are  interesting.  Problems  will  defined  on  some  particular  classes  in  Sections  4.1 
and  5.  A  class  can  be  identified  using  a  5-tuple  in  which  each  component  specifies  one  of  the  following  details 
about  the  networks  in  the  class. 

1.  Restrictions  on  parallel  arcs:  Parallel  arcs  may  or  may  not  be  allowed  in  the  tree  part  of  the  network. 
This  restriction  refers  to  the  directed  arcs  in  G,  not  to  the  undirected  edges  mentioned  in  the  definition 
of  an  STT  network.  Excluding  parallel  arcs  does  not  preclude  the  use  of  anti-parallel  arcs.  This 
restriction  does  not  apply  to  arcs  to  or  from  node  0. 

2.  Restrictions  on  arc  flow  capacities:  The  arcs  of  a  network  can  have  either  unit  capacities  or  general 
non-negative  integer  capacities. 

3.  Restrictions  on  arc  flow  multipliers:  Multipliers  can  either  be  unity,  as  in  an  ordinary  network,  or 
general  non-negative  integers,  as  in  a  generalized  network. 

4.  Restrictions  on  node  demands:  A  network  may  consist  of  all  transshipment  nodes,  all  demand  nodes, 
all  supply  nodes,  or  a  mixture.  These  descriptions  do  not  apply  to  node  0. 

5.  Restrictions  on  node  excesses  for  feasible  flows:  A  flow  may  be  required  to  be  sufficient,  deficient,  or 
exact.  Note  that  this  is  a  condition  on  acceptable  solutions,  whereas  the  first  four  conditions  refer  to 
characteristics  of  the  graph  structure. 

The  possible  values  for  each  component  of  the  5-tuple  are  listed  in  Figure  1.  For  example,  the  class  of 
networks  C  —  (Al/CG/MG/N M/ FS)  consists  of  STT  networks  with  no  directed  parallel  arcs,  capacities 
and  multipliers  that  can  be  any  non-negative  integers,  both  supply  and  demand  nodes,  and  in  which  the 
flow  to  each  node  must  be  sufficient.  This  class  of  STT  networks  is  discussed  further  in  Section  4.1. 

A  flow  X  G  is  feasible  if  it  satisfies  the  flow  requirements  (the  fifth  component  of  the  STT  class 

definition)  and  if,  in  addition,  x  <  u.  The  set  of  feasible  flows  for  an  instance  constitutes  the  feasible  set 
S  for  that  instance;  the  integer  program  that  specifies  the  set  is  generally  not  written  out  explicitly.  The 
collection  of  feasible  sets  S  for  the  networks  in  a  class  of  STT  networks  is  the  family  of  feasible  sets  ^  for 
an  STT  Network  Flow  problem  defined  on  the  class. 

This  section  has  described  feasible  sets  for  STT  Network  Flow  problems.  In  order  to  specify  an  STT  Network 
Flow  optimization  or  feasibility  problem  completely,  other  parts  of  the  problem  must  be  specified  as  well:  a 
family  T  of  r-criteria  objective  functions  defined  on  the  vector  x  of  flows,  an  r-dimensional  direction  vector 
w,  and,  for  a  feasibility  problem,  an  r-dimensional  target  vector.  The  next  section  discusses  some  appropriate 
classes  of  objective  functions,  and  Section  4.1  presents  a  VPP  algorithm  for  a  problem  defined  on  a  particular 
class  of  STT  networks. 
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Arcs: 

A1 

AP 

No  parallel  arcs  are  allowed  in  the  tree  part  of  the  network. 
Parallel  arcs  are  allowed. 

Capacities: 

Cl 

Unit  flow  capacities  (w  =  1). 

CG 

General  flow  capacities  in  Z'^. 

Multipliers: 

Ml 

Unit  flow  multipliers  (fi  =  1). 

MG 

General  flow  multipliers  in  Z'^ . 

Nodes: 

NO 

Network  may  contain  only  transshipment  nodes  (6  =  0). 

ND 

Network  may  contain  only  demand  nodes  (6  >  0). 

NS 

Network  may  contain  only  supply  nodes  (6  <  0). 

NM 

Network  may  contain  a  mixed  node  set  (no  restriction  on  6). 

Flows: 

FS 

Flow  into  each  node  must  be  sufficient. 

FD 

Flow  into  each  node  must  be  deficient. 

FE 

Flow  into  each  node  must  be  exact. 

FU 

Flow  excesses  are  unrestricted. 

Figure  1:  Components  of  STT  network  class  specification. 

3  Objective  Functions 

The  problems  discussed  in  the  remainder  of  this  paper  are  distinguished  from  each  other  primarily  by  the 
structures  of  their  feasible  regions.  Before  looking  at  these  feasible  regions,  however,  some  attention  should 
be  paid  to  the  objective  functions  that  are  used.  This  section  examines  several  kinds  of  objective  functions 
that  arise  frequently  in  practice. 

This  section  does  not  deal  with  an  important  aspect  of  objective  functions:  whether  they  are  to  be  minimized 
or  maximized.  The  difference  is  crucial  when  considering  approximation  algorithms.  Consider  the  Traveling 
Salesman  problem  on  complete  graphs  without  the  triangle  inequality.  The  minimum  cost  version  has  no 
polynomial  time  £-approximate  algorithm  for  any  e  >  0  unless  V  =  A/’'P[SG76].  A  tour  of  length  at  least 
one-half  the  optimum  can  easily  be  found  in  the  maximization  case,  however,  by  extending  a  maximum 
weight  matching. 

Another  example  is  the  Multi-Dimensional  Binary  Knapsack  problem  discussed  in  [S094].  As  mentioned 
there,  the  maximization  form  has  no  FAS  unless  V  =  NV .  The  minimization  form,  however,  can  be  solved 
exactly  and  quickly  by  setting  a;  =  0,  that  is,  leaving  the  knapsack  empty. 

Section  3.1  generalizes  the  notion  of  monotonically  increasing  functions.  This  generalization  is  applied  to 
additively  separable  functions  in  Section  3.2.  VPP  reductions  between  problems  with  difference  function 
classes  are  derived  in  Section  3.3. 


3.1  Monotonicity 

Many  common  objective  functions  are  monotonically  increasing  or  monotonically  decreasing.  A  function 
/  :  ^  Z’’"*'  is  monotonically  increasing  if  for  ar,  j/  E  2"'*',  x  <y  implies  that  f  {x)  <  f  {y).  The  “  <  ” 
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relation  is  applied  on  a  component-by-component  basis.  The  function  /  is  monotonicaJly  decreasing  if  the 
function  g  =  — /  is  monotonically  increasing. 

Results  for  these  kinds  of  functions  have  been  derived  in  the  single-criterion  case.  Restricting  attention  to, 
say,  monotonically  increasing  objective  functions  in  the  multi-criteria  case,  however,  is  not  as  useful,  because 
different  components  of  an  objective  function  can  be  optimized  in  different  directions.  One  way  to  generalize 
single-criterion  monotonic  functions  is  to  require  that  the  values  of  all  components  of  the  objective  function 
simultaneously  improve  or  get  worse  when  the  components  of  the  argument  increase,  as  illustrated  by  the 
following  definitions. 

Definition  1  An  objective  function  /  :  is  order-preserving  w.r.t.  a  direction  vector  w  G  fi'’  if 

either  of  the  following  two  conditions  holds: 


'  Wi  =  “  <  ” 

/»•  (»)  <  fi  {y) 

(a) 

V  a;,  j/  G  ,  x  <y  for  i  =  1 , . . 

.,r,  < 

Wi  =  “  >  ” 

fi  (®)  >  fi  (y) 

,  Wi  =  “  =  ” 

fi  (a;)  <  fi  (y) 

Ui  —  5. 

fi  {x)  <  fi  {y) 

(b) 

'i  x,y  S""*",  X  <y  =>  for  «  =  1, . . 

.,r,  ^ 

Wj  —  ,> 

=» 

fi  («)  >  fi  {y) 

,  s.  —  “  —  ” 

t  Wi  _  - 

fi  {x)  >  fi  {y) 

The  point  is  that  the  value  of  each  component  of  an  order-preserving  function  changes  in  accordance  with  the 
corresponding  component  of  the  direction  vector.  Component  functions  corresponding  to  direction  vector 
entries  of  “  =  ”  can  behave  as  either  “  <  ”  or  “  >  ” ,  but  must  be  consistent  in  this  regard  over  the  entire 
feasible  set.  Furthermore,  all  component  functions  that  correspond  to  “  =  ”  entries  must  behave  the  same 
way. 

The  definitions  of  “monotonically  increasing”  and  “order  preserving”  could  have  been  stated  for  functions 
that  map  from  a  subset  S  C  2”"^  rather  than  from  the  entire  space  Z’^'^.  Doing  so  would,  however,  have 
required  a  large  overhead  in  notation  to  accommodate  the  function  classes  defined  later  in  this  section.  The 
more  complicated  formulation  would  not  add  to  the  insight  provided  by  the  results  in  this  paper,  nor  would 
it  extend  the  applicability  of  the  results  in  practice.  Because  of  these  considerations,  the  present  definitions 
are  used  instead  of  more  general  definitions. 

Definition  2  An  objective  function  /  is  order-reversing  w.r.t.  a  direction  vector  w  if  the  function  g  —  —f  is 
order-preserving  w.r.t.  w,  and  /  is  order-monotone  w.r.t.  a  direction  vector  w  if  /  is  either  order-preserving 
or  order-reversing  w.r.t.  to.  ■ 


-5- 


Multi-Criteria  Flow,  Knapsack,  and  Scheduling  Problems 


Definition  3  The  following  classes  of  objective  functions  are  defined  for  a  direction  vector  w: 

V  (w)  =  {  /  :  /  is  order-preserving  w.r.t.  w  }  , 

'R.{u!)  =  {  /  :  /  is  order-reversing  w.r.t.  w  }  , 

At  (w)  =  V  (u>)  LI  IZ  (u)  . 


The  notions  of  closure  under  scaling  and  under  box  constraints  are  introduced  in  [S094].  For  any  direction 
vector  u,  each  of  the  sets  V  (w),  TZ  (w),  and  A1  (w)  is  closed  under  scaling  and  under  box  constraints.  The 
following  lemma  describes  the  relationship  between  the  sets  V  (w)  and  TZ{ui). 

Lemma  1  Let  S  be  an  n-dimensional  feasible  set  containing  a  “largest"  element  u  that  satisfies  the  condition 
®  €  5"  =>  (u  —  x)  £  S.  Let  f  :  S  -2’’+  be  a  function  and  define  the  function  f  :  S  ^  2'’+  by 
f{x)  =  f  {u  —  x).  Then  f  EV  (w)  if  and  only  if  f  E  TZlu)). 

Furthermore,  a  component  function  fi  that  corresponds  to  LOi  =  “  —  behaves  as  though  w,-  =  “<  «  ^ 

according  to  pari  ( a)  in  the  definition  of  order-preserving)  if  and  only  if  the  function  f-  behaves  as  though 

Ui  —  “>  ”  (i.e.,  according  to  part  (b)  in  the  definition  of  order-preserving). 

Proof.  Suppose  that  /  G  T{co).  Let  x,y  E  S  he  given,  with  x  <  y;  then  (u  -  a;),(u  —  y)  E  S,  and 

{u  —  x)  >  (u  —  y).  For  each  i  =  1, . . . ,  r,  consider  the  possible  values  of  Wj.  If  w,-  ==  “  <  ” ,  then 

//(«)  =  fi{u-x)>  fi  {u-y)  =  fl(y)  . 


Analogous  statements  are  true  for  Wj  G  Therefore  f  ElZ  (w). 

The  proof  of  the  converse  statement  is  similar.  ■ 


3.2  Separability 

The  necessary  and  sufficient  conditions  in  [S094]  for  the  existence  of  an  FAS  apply  to  problems  with  quite 
general  objective  functions.  Often,  though,  VPP  algorithms  use  dynamic  programming,  which  limits  the 
kinds  of  objective  functions  that  can  be  accommodated  easily.  In  order  to  be  amenable  to  the  application 
of  dynamic  programming,  a  problem  must  typically  have  a  separable  objective  function.  The  most  common 
separable  objective  functions  are  additively  separable,  that  is,  they  can  be  written  as  f  [x)  —  fj  {xj). 

Functions  of  the  form  f  (x)  -  n"=i  fj  (xj)  and  f  (x)  =  min{/j  (xj) ,  j  =  1, . . .  ,n}  are  also  used[Whi69, 
section  3.2].  This  section  introduces  additively  separable  functions  that  are  order-monotone. 

The  reason  that  separable  objective  functions  are  common  in  dynamic  programming  applications  is  that 
the  value  of  the  function  can  be  determined  incrementally  by  optimizing  the  value  of  one  variable  at  a 
time.  Although  the  focus  here  is  additively  separable  functions,  fairly  general  definitions  of  separability  are 
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presented  in  [YS81,  VK82,  Har83].  The  problems  that  are  amenable  to  solution  by  dynamic  programming 
are  characterized  in  [Bel61],  as  cited  in  [LH90],  as  those  which  can  be  formulated  as  Markov  processes.  Note 
that  the  conditions  of  monotonicity  and  nonbounded  dominance  in  [YS81]  are  satisfied  by  the  objective 
functions  discussed  in  [S094]. 

Some  attempts  have  been  made  to  extend  dynamic  programming  to  problems  with  objective  functions 
that  are  not  separable.  One  approach  is  to  consider  history-remembering  policies[Kre77,  Hen85].  Another 
approach  is  to  approximate  the  non-separable  functions  with  separable  functions[Sni87].  These  approaches 
are  not  pursued  in  this  work. 

Definition  4  A  function  f  :  S  ^  ,  where  S  C  is  additively  separable  if  there  exist  functions  /i, 

. . . ,  fn  such  that  for  all  x  £  S,  f{x)  =  X3"=i  fj  (®j)-  ■ 

The  discussion  in  Section  3.1  can  be  specialized  to  this  kind  of  function. 

Definition  5  For  a  direction  vector  ca  6  O'", 

AV  (w)  =  {  /  G  P  (w)  :  /  is  additively  separable  }  , 

ATI  (m)  =  {  /  G  77  (w)  :  /  is  additively  separable  }  , 

AM{io)  =  AV{w)uA1l{w)  . 


The  notions  of  quasi-closure  under  scaling  and  under  box  constraints,  as  well  as  of  scaling  neighbor  and  box 
constraint  neighbor,  were  introduced  in  [S094]  along  with  the  notions  of  closure  mentioned  in  Section  3.1. 
None  of  the  sets  in  Definition  5  is  closed  under  scaling  or  under  box  constraints,  but  for  any  direction  vector 
u,  each  is  quasi-closed  under  scaling  and  under  box  constraints. 

Scaling  and  box  constraint  neighbors  of  a  function  /  G  AAd  (w)  are  easy  to  find.  Write  /  as  /  =  fj  (®l)i 
where  fj  is  an  r-component  function.  For  each  j  —  1, . . . ,  n  and  i  —  1, . . .  ,r,  let  fij  be  the  i  component 
function  of  fj.  Then  one  scaling  neighbor  of  /  w.r.t.  v  is  the  function  g  with  r  component  functions 
9i  —  Sj=i  /'^*J !  *  =  1>  •  •  • !  J*-  One  box  constraint  neighbor  of  /  w.r.t.  M  is  the  function  h  with  r 

component  functions  hi  —  {fi,i  iZi)  ^  M,},  i  =  1, . . . ,  r. 

The  functions  that  constitute  the  terms  of  an  additively  separable  function  behave  similarly  to  the  entire 
function  with  regard  to  order  monotonicity,  as  shown  by  the  following  lemma. 

Lemma  2  Let  S  C  be  closed  under  the  “  <  ”  relation.  That  is,  for  any  x,y  £  Z'^^ ,  with  x  <  y,  if 
y  £  S  then  x  £  S.  Suppose  that  a  function  f  :  S  Z^'^  is  defined  by  f  (x)  =  fj  i^j)-  Then 

(a)  f  £  AV  (w)  if  and  only  if  for  j  =  1, . . .  ,n,  fj  £V  (w),  with  the  same  behavior  on  =  “=  ”  as  f,  and 

(b)  f  £  ATi{w)  if  and  only  if  for  j  =  1,  ■  ■  ■  ,n,  fj  £  77  (w),  with  the  same  behavior  on  Wi  —  “=  ”  as  f. 


-7- 


Multi-Criteria  Flow,  Knapsack,  and  Scheduling  Problems 


Proof.  Only  the  proof  for  (a)  is  given;  the  proof  for  (b)  is  similar. 

Necessity:  Suppose  that  /  €  AV  (w).  Choose  x,y  £  S  with  x  <  y,  and  focus  on  some  criterion  i.  Suppose 
that  u)i  =  “  <  ” ;  the  proof  is  similar  for  Wi  . 

Since  /  £  AV  (w),  fi  (a:)  <  /j  (y).  Suppose  that  fj  does  not  behave  as  claimed  for  some  j  £  {I,...,  n].  Let 
X  =  (0, . . . ,  0,  Kj,  0, . . . ,  0)  and  y  =  (0, . . . ,  0,  %■,  0, . . . ,  0),  where  both  x  and  y  are  n-vectors  with  the  non-zero 
elements  in  component  j. 

By  assumption,  x,y  £  S.  But  f  (x)  =  fj  {x)  and  /(y)  =  /_,■  (y),  contradicting  the  assumption  that  does 
not  behave  the  in  the  same  way  as  /. 

Sufficiency:  Suppose  that  for  j  =  1, . . . ,  n,  fj  £  V  (w),  all  with  the  same  behavior  on  uji  =  “  =  and  that 
x,y  £  S  with  X  <  y.  For  any  i  =  suppose  that  Wj  =  “  <  Then  fi{x)  =  I3"=i /bi  (®i)  < 

Z)"=i  fij  iV})  =  fi  {y)-  Analogous  statements  are  true  for  w,-  So  f  £  AV  (w),  with  the 

same  behavior  on  wj  =  “  =  ”  as  the  functions  fj .  ■ 

3.3  Reductions  for  Classes  of  Objective  Functions 

Reductions  can  be  used  to  show  relationships  between  classes  of  objective  functions.  The  objective  function 
classes  defined  in  Section  3  are  of  particular  interest  in  the  present  work. 

Consider  first  the  classes  of  order-preserving  and  order-reversing  objective  functions.  Problems  defined  with 
these  function  classes  are  similar  in  some  respects.  This  is  particularly  true  when  the  direction  vectors  are 
in  that  is,  the  set  of  direction  vectors  {ui^’=  :  r  =  1, . . .}.  The  relationship  between  'P(w’’>=), 
and  M  for  any  r,  is  described  in  the  following  theorem. 

Theorem  1  For  any  family  of  feasible  sets,  the  following  problems  are  VPP  equivalent  for  any  k: 

1.  (^,  Ad  (w”’") ,  w’’’=,  Teas)  ,  r  =  l,2,...,k. 

2.  {^,V{oj^>=),u^'=,Feas),  r  =  l,2,...,k. 

3.  ,  Feas) ,  r  =  l,2,...,k. 

Proof.  For  any  r,  both  (^,P  (w’''=) ,  Feas)  and  Feas)  are  the  same  as,  and  so  are 

VPP  equivalent  to,  (■$,  Ad  {w’’’-)  Feas). 

Assume  for  the  moment  that  r-criteria  problems  can  be  VPP  reduced  to  single-criterion  problems,  i.e.,  that 
Ad  ,  Feas)  ocvpp  (^,  Ad  Feas);  this  assertion  will  be  proved  below.  The  problem 

(^,Ad  (w^’~)  ,a;^>“',Feas)  can  be  VPP  reduced  to  each  of  the  other  listed  problems  as  follows: 

1.  An  r-criteria  function  can  be  considered  as  an  (r  +  l)-criteria  function  with  an  (r  -f  1)®*  criterion 
that  is  identically  zero,  so  Ad  (w'’’^)  C  Ad  (a)’'+b=).  Use  Lemma  1  in  [S094]  to  conclude  that 
(^,  Ad  (wb=)  Feas)  oCvpp  (^,  Ad  ,  w’’’-,  Feas). 
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2.  As  noted  above,  (^,Ad  ,a;b-,Feas)  oCvpp  {^,1^  ,wb-^Feas).  The  further  VPP  reduction 

to  {^,V  ,0;'’’=,  Feas)  follows  as  in  (1). 

3.  The  reduction  to  (^,  7^  Feas)  is  similar  to  that  in  (2). 

The  theorem  is  proved  by  establishing  the  earlier  assertion: 

,0;'’’=,  Feas)  oCvpp  (^,  Ad  (wb=)  ,w^>=,Feas): 

Let  I  =  {S,  ,M)  e  (^,  Ad  ,a;’'>=,Feas)  be  given.  Define  the  function  h  •.  S  —*  Z'^  by  h{x)  = 

Z)i=i  /«■  (®)  (/)  -b  l]*.  The  function  h  can  be  computed  in  VPP  time.  For  any  v  G  Z''+,  f  {x)  =  u  if 

and  only  if  h{x)  =  Wj  [My  (/)  -f  l]  .  Furthermore,  since  each  direction  component  iJl‘~  has  the  value 
“  =  ” ,  either  each  component  function  fi  is  monotonically  increasing  or  each  is  monotonically  decreasing; 
hence  so  is  /i.  This  means  that  h  G  Ad  (wb=)_  Therefore  the  instance 
(^,Ad  (wb=)  jwi>=,Feas).  ■ 

Other  relationships  become  apparent  when  the  directions  of  optimization  are  reversed. 

Definition  6  The  dual  of  a  direction  vector  w  is  obtained  by  changing  “  <  ”  entries  in  w  to  “  >  ”  and 
vice  versa.  ■ 

Example  1  (Dual  direction  vector)  Here  is  an  example  of  taking  the  dual  of  a  direction  vector: 


/  <  ^ 

3LJ 

i  -  ^ 

> 

= 

< 

\  = ) 

V  =  / 

Also  note  that  for  any  r,  =  u’'’~ .  ■ 

Lemma  3  For  any  family  of  feasible  sets  and  direction  vector  w, 

(a)  (^,7’(w),w,Feas)  =vpp  (^,  7^  (w)  ,  w®  ,  Peas)  . 

(h)  (^,'P  (w)  ,w^.  Peas)  =vpp  (^,72.(w)  ,w,  Peas)  . 

Proof.  Only  the  “only  if”  part  of  (a)  is  proved;  the  “if”  part  of  (a)  and  both  parts  of  (b)  can  be  proved 
analogously.  Let  (S', /,  w,  M)  G  (^,  P  (w)  ,w,  Feas)  be  given.  Let/'=min{/,  M-fer}.  SinceP(w)  is  closed 
under  box  constraints,  f  EV  (w),  and  (/'  (x) ,  M)  if  and  only  if  (/  (a?) ,  M). 

Define  g  —  M  +  Cr  —  f;  then  g  G  P  (w).  For  any  x  E  S,  (/'  (x) ,  M)  if  and  only  if  D^JD  {g  (a;) ,  e^). 
So  solving  the  instance  {S,g,u^,er)  G  (^,P  (w) ,  Feas)  yields  a  solution  to  the  instance  (S, /,w,M)  G 
(^,  P  (w) ,  w,  Feas).  ■ 

Similar  relationships  are  true  when  consideration  is  restricted  to  additively  separable  functions.  The  proofs 
of  the  following  theorem  and  lemma  are  analogous  to  the  proofs  of  Theorem  1  and  Lemma  3. 


-  9  - 


Multi-Criteria  Flow,  Knapsack,  and  Scheduling  Problems 


Theorem  2  For  any  family  ^  of  feasible  sets,  all  the  following  problems  are  VPP  equivalent: 

1.  ,w'^’~,Feas) ,  r=l,2, ...  . 

2.  {■9,AViw’-’=),u>^’=,Feas),  r=l,2,.... 

3.  ('9,Aniu''’=),w^’=,Feas),  r=l,2,.... 

Lemma  4  For  any  family  W  of  feasible  sets  and  direction  vector  u>  €  Q*", 

(a)  (^,AP{u>),w,Feas)  =vpp  {^9 ,  ATI  (ao) ,  uj^  ,  Feas)  . 

(b)  ifS ,AV  {u)  ,Feas)  =vpp  ,A'll{u))  ,uj,  Feas)  . 

4  Some  Simple  STT  Network  Flow  Problems 

The  class  of  STT  networks  C  =  (Al/CG/MG/NM/ FS)  is  described  in  Section  2.  Section  4.1  focuses  on  a 
canonical  subset  C'  of  that  class.  A  VPP  algorithm  is  presented  for  a  feasibility  problem  defined  on  C'  and 
is  used  to  prove  the  existence  of  an  FAS  for  the  optimization  form  of  that  problem.  VPP  reduction  is  used 
to  explore  the  relationship  between  C  and  C'  in  Section  4.2. 

4.1  A  VPP  Algorithm 

A  network  in  C  is  in  canonical  form  if  a  single  arc  connects  node  0  to  each  other  node,  each  other  connected 
pair  of  nodes  is  connected  by  a  single  pair  of  anti-parallel  arcs,  and  the  tree  part  of  the  network  is  connected. 
More  formally,  a  network  G  G  C  is  in  canonical  form  if  it  satisfies  the  following  two  conditions: 

1.  Vv  G  V  —  {0},  A  contains  a  single  arc  (0,  v)  and  no  arc  (v,  0). 

2.  The  tree  part  of  G  is  connected.  If  nodes  v  and  w  are  adjacent,  then  A  contains  a  single  arc  (v,  w)  and 
a  single  arc  {w, «). 

The  collection  of  networks  of  C  in  canonical  form  is  denoted  C . 

A  tree  T,  distinct  from  but  related  to  the  tree  part  of  G,  is  eissociated  with  G  for  purposes  of  describing 
the  VPP  algorithm.  The  tree  T  has  nodes  {1, . .  .,m};  node  v  in  T  corresponds  to  node  v  in  G.  For  any 
v,w  E  {1, . . . ,  m},  nodes  v  and  wofT  are  connected  by  an  arc  if  and  only  if  nodes  v  and  w  are  adjacent  in 
G.  The  directions  of  the  tree  arcs  are  specified  next. 

In  a  postorder  traversal  of  a  tree,  each  subtree  of  the  root  is  traversed  (recursively,  in  postorder)  and  then 
the  root  is  visited [Knu73,  RND77].  Choose  an  arbitrary  node  p  as  the  root  of  T.  Label  the  nodes  according 
to  any  postorder  ordering  w.r.t.  the  root  p,  representing  the  label  of  node  v  by  /  (v).  In  particular,  I  {p)  —  n. 
The  VPP  dynamic  programming  algorithm  to  be  described  processes  the  nodes  of  N  -  {0}  in  increasing 
label  order.  A  good  overview  of  this  approach  can  be  found  in  [JN83],  from  which  some  notation  is  adapted. 
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Direct  all  the  arcs  of  T  away  from  p  and  let  be  the  number  of  children  of  node  v  in  T.  This  number  is  also 
called  the  outdegree  of  v.  For  each  node  v  and  for  j  =  1, . .  .,dv,  let  T[v,  j]  be  the  subtree  of  T  defined  by 
V,  the  first  j  children  of  v,  and  all  the  descendants  of  those  children,  where  the  children  are  ordered  by  label 
value.  For  example,  T[i;,0]  is  just  the  node  v  and  T[p,dp]  —  T.  Let  G{v,i]  be  the  subgraph  of  G  defined  by 
node  0  and  the  nodes  of  G  that  correspond  to  the  nodes  of  T[v,  j]. 

A  VPP  algorithm  that  solves  the  feasibility  problem  (C'jAV  ,w^'-,Feas)  will  be  presented.  For 

purposes  of  stating  the  algorithm,  say  that  a  flow  a;  in  G  is  feasible  w.r.t.  G[v,]\  if  it  satisfies  the  following 
three  conditions: 

1.  The  flow  into  each  node  of  G[v,j],  except  possibly  into  node  v  itself,  is  sufficient. 

2.  If  {w,w')  ^  Glv,j],  then  =  0. 

Z.  X  £  and  x  <  u. 

Define  9(v,j,k)  to  be  the  maximum  excess  flow  into  node  v  over  all  flows  x  that  both  are  feasible  w.r.t. 
G\v,j\  and  have  f  {x)  =  k.  The  algorithm  attacks  the  flow  problem  indirectly  in  that  the  computational 
focus  is  the  maximum  node  excesses  9  rather  than  the  objective  function  /.  This  approach  works  because 
of  the  following  lemma. 

The  idea  behind  this  lemma  is  to  examine  9{p,dp,  M).  The  flows  considered  in  determining  the  value  of 
9(p,  dp,  M)  are,  by  the  definition  of  9,  feasible  on  the  subtrees  of  the  root  node  p.  Such  a  flow  is  feasible  for 
the  entire  graph  G  if  and  only  if  the  excess  flow  into  p  is  sufficient,  that  is,  if  and  only  if  9{p,  dp,  M)  >  0. 

Lemma  5  An  instance  (G,  /,  “=  ”,  M)  of  the  problem  {C'  ,A'P  Feas)  has  a  solution  if  and  only 

if9(p,  dp,M)  >  0. 

Proof. 

Sufficiency:  Suppose  that  the  instance  (G,  /,  “  =  "  ,M)  has  a  solution  x.  Then  a;  is  a  feasible  flow,  so 

1.  The  flow  into  each  node  of  G{p,  dp]  is  sufficient  (including  the  flow  into  p). 

2.  All  arcs  are  in  G[/9,  df]. 

Z.  X  £  and  x  <u. 

So  X  is  feasible  w.r.t.  G[p,df\  and  f  (x)  =  M.  Furthermore,  since  the  flow  into  p  is  sufficient,  Cp  (x)  >  0. 
Since  9{p,  dp,  M)  is  the  maximum  of  a  set  that  contains  Cp  (x),  it  must  be  that  9{p,  dp,  M)  >  0. 

Necessity:  Suppose  that  9(p,dp,  M)  >  0.  This  implies  the  existence  of  a  flow  x  that  is  feasible  w.r.t. 
Cr\p,  df\  =  G  and  satisfies  /  (a:)  =  M .  Furthermore,  Cp  (a;)  =  9{p,  dp,  M)  >  0,  so  the  flow  into  every  node  is 
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for  i—  1 , . . . ,  n 

Let  V  be  the  node  with  /  (v)  =  i. 
for  j  =  0, . . . ,  dv 

for  fe  =  0, . . . ,  M 

(*)  Compute  d(n,j,fe). 

if  $(p,  dp,  M)  >  0 

then  X  is  the  desired  solution, 
else  the  problem  is  infeasible. 


Figure  2:  The  algorithm  for  (C',.4'P  ,w^’~,Feas). 

sufficient.  Therefore  a;  is  a  feasible  flow  that  solves  the  instance  (G,  /,  “  =  M).  ■ 

The  algorithm  for  solving  (^CjAV  ,w^’“,Feas)  appears  in  Figure  2.  A  detailed  description  of  the 

main  computation  step,  labeled  (*),  appears  next.  This  is  followed  by  a  proof  that  the  algorithm  is  VPP  for 
{C',AV  ,w^'=,Feas). 

The  computation  for  (=t=)  has  two  cases,  depending  on  whether  j  =  0  or  j  >  0.  When  j  =  0,  perform  the 
computation  by  solving  the  following  program. 

6{v,  0,  k)  =  max  po,vXo,v  —  bv 

2/  0 ,  v 

s.t.  /o,v  (^o,v)  ~  k 
^0,v  ^  ^0,v 

xo.v  e  2+ 


When  j  >  0,  let  w  be  the  child  of  v,  and  perform  the  computation  by  solving  the  following  program. 


max  0{v,j  —  l,a)  +  fi,^^yXn,^v-Xv,w 

Ol  f  ^  y  'y  y  By  ^  W  yV 

(2a) 

S.t.  0(^Wj  dyj  J  yW^V  ,W  ^VJyV  ^ 

(2b) 

fVyW  ~  7 

(2c) 

fWyV  i^WyV)  —  ^ 

(2d) 

(2e) 

^VyiV  ^  '^Vy'HU 

(2f) 

^WyV  ^  '^W  ,V 

(2g) 

/^J  T)  ^VyW:  ^iVyV  € 

(2h) 

Theorem  3  The  algorithm  in  Figure  2  is  VPP  for  the  feasibility  problem  {C' ,A'P  Feas). 

Proof.  The  proof  has  two  parts:  the  algorithm’s  correctness  and  its  running  time. 

Correctness:  In  the  CEise  that  j  =  0,  the  maximum  excess  on  G[v,  0]  must  be  computed.  This  is  just  the 
subgraph  of  G  induced  by  the  nodes  0  and  v,  which,  because  G  G  C' ,  contains  a  single  arc,  (0,i;).  The 
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quantity  iio,v^o,v  —  is  the  excess  flow  into  node  v  on  G[v,  0]  when  the  flow  on  arc  (0,t;)  is  The 

program  specified  by  Equations  1  maximizes  this  excess  subject  to  the  feaisibility  conditions,  and  so  correctly 
computes  6{v,0,k). 

In  the  case  that  j  >  0,  the  maximum  excess  on  G[v,  j]  must  be  computed.  Note  that  the  nodes  are  considered 
in  increasing  label  order,  the  postorder  label  of  a  node  is  greater  than  the  labels  of  its  children,  and  the 
children  of  a  node  are  processed  in  increasing  label  order;  therefore  the  values  of  that  appear  in 

Equations  2a  and  2b  are  computed  before  they  are  used. 

The  structure  of  Equations  2  comes  from  decomposing  the  flow  in  G[v,  j]  and  the  value  k  into  four  parts:  the 
flow  in  G[v,j  —  1],  with  value  a;  the  flow  in  G[w,  dto],  with  value  ,0;  and  the  flows  on  arcs  [v,  w)  and  (w,  v), 
with  values  7  and  S.  The  total  required  value  of  k  can  be  partitioned  in  this  way  because  additively  separable 
objective  functions  are  used. 

The  objective  function,  Equation  2a,  expresses  the  value  of  the  excess  flow  into  node  v,  given  the  flows  on 
arcs  (?;,  w)  and  (w,  v)  and  in  the  first  (j  —  1)  subtrees  of  node  v,  and  subject  to  the  flow  in  G[v,  j  —  1]  having 
value  a.  This  is  maximized  to  determine  9{v,  j,  k). 

Equation  2b  ensures  that  the  flow  into  node  w  is  sufficient;  recall  that  it  need  not  have  been  sufficient  in 
G[w,dw\,  and  anyhow  the  arcs  (u,  tn)  and  {w,v)  are  not  in  G[w,dw].  The  value  contributed  by  the  flow  in 
G[w,  d^]  is  constrained  to  be  /?.  Since  Equation  2b  expresses  the  feasibility  of  the  flow,  the  excess  might  as 
well  be  as  large  as  possible,  so  0(w,dw,^)  is  assumed. 

Equations  2c  and  2d  constrain  the  flow  values  on  arcs  {v,w)  and  (m, «).  Equation  2e  ensures  that  the  four 
parts  of  the  flow  add  up  to  fc,  and  Equations  2f  through  2h  are  needed  to  guarantee  feasibility. 

So  the  program  specified  by  Equations  2  determines  9{v,j,k).  The  final  value  computed  by  the  algorithm 
is  9{p,  dp,  M).  If  this  value  is  non-negative,  then  the  necessity  part  of  the  proof  of  Lemma  5  shows  that  the 
value  of  X  found  by  the  algorithm  is  a  solution  to  the  instance  {G,  /,  “  =  ” ,  M).  Similarly,  if  9{p,  dp,  M)  <  0, 
then  Lemma  5  implies  that  the  instance  has  no  solution.  Therefore  the  algorithm  is  correct. 

Running  time:  The  computation  step  (+)  is  executed  only  M  ■  \A\  times.  The  rest  of  this  proof  shows  that 
each  execution  can  be  solved  in  VPP  time.  This  yields  the  conclusion  that  the  algorithm  runs  in  VPP  time. 

A  call  with  j  —  0  requires  solution  of  the  program  expressed  by  Equations  1.  For  a  node  v,  this  can  be 
solved  with  O  (log  +  1))  evaluations  of  /o,„  (•)  by  using  binary  search  over  possible  values  of  *o,k  ,  because 
fo,v  G  and  fo,v  is  therefore  either  monotonically  increasing  or  monotonically  decreasing. 

A  call  with  j  >  0  requires  solution  of  the  program  expressed  by  Equations  2.  Only  O  (fc^)  distinct  4-tuples 
(a,  /?,  7,  S)  satisfy  Equations  2e  and  2h,  so  if  the  program  can  be  solved  for  each  4-tuple  in  VPP  time,  then 
the  entire  algorithm  will  run  within  the  allowed  time. 
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Fix  values  of  a,  jS,  j,  and  S.  The  value  of  Xv,w  must  satisfy 

fv,w  —  T 

,W  <  Uv  yW 

Xv,w  €  2'^ 


(3) 


Let  and  be  the  smallest  and  largest  values  that  satisfy  Equations  3;  they  can  be  found  with 
0(log(«„_„  +  1))  evaluations  of  (•)  by  binary  search,  because  6  AV  Similarly,  the  corre¬ 

sponding  values  of  and  Cw,v  for  can  be  found  with  0(log(w„,  „  -|-  1))  evaluations  of  (•).  The 
following  program  yields  the  optimal  values  of  and  given  (a,/?,  7, (5). 


^  V  .W  ]  w  ,v 


f^W 

,v  " 

,w 

,w 

VI 

3 

^v. 

,w  ^  ^v,w 

Ay; 

VI 

XjjD 

,V  ^  Cw  jV 

J 

W  )  ^WjV 

(4) 


To  see  that  this  program  can  be  solved  in  VPP  time,  consider  the  following  cases: 

1-  f^V  ,W  —  0*  Set  ,W  •  If  ^  f^VjW^VfW  -t-  9{'w,dyj,P)  —  by,,  then  the  program  is  infeasible. 

Otherwise,  set 

Xyj  jU  —  min  f  P'V  ,W  Xy  ^y,  -f  0{w ,  dyj  ,  by,  }  . 

2-  fJ.w,v  —  0,  fiy^w  >  1:  Set  Xy,^y  =  If  \bw  —  d{w,dy,,  Xy,^y]  /  fiy^y,  >  ^y^y,,  then  the  program  is 

infeasible.  Otherwise,  set 


Xy^yj  -  max  ■{  ,  [6^,  9^W  ,  dy,  ,  -j-  j  j fly^y,  ^ 


3.  fiy^y,,  Hw^y  >  1\  Supposo  that,  in  some  feasible  solution  of  Equations  4,  both  Xy^y,  <  ^y^y,  and  Xy,^y  < 
^w^y.  Incrementing  each  variable  by  one  would  yield  a  new  feasible  solution  with  objective  value  at 
least  as  good  as  the  value  of  the  previous  solution.  Therefore,  attention  can  be  restricted  to  solutions 
in  which  at  least  one  of  these  variables  is  at  its  upper  bound.  Each  case  is  considered,  and  the  one 
with  the  larger  objective  function  value  is  used. 

(aj  Xy^y,  —  ^y  ,y}  .  Set  Xyj^y  as  in  case  1. 

^bj  Xy/^y  —  ,v  ‘  Set  Xy^yj  ss  in  case  2. 


Apply  Lemma  4  and  Lemma  3  in  [S094]  to  Theorem  3  in  order  to  conclude  that  the  optimization  problem 
corresponding  to  (C'jAV  (wb=)  ,w^>=,Feas)  has  an  FAS. 
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Corollary  1  The  optimization  problem  {C' ,A'P  Opi)  has  an  FAS. 

4.2  Using  Reduction 

The  previous  section  presented  an  explicit  VPP  algorithm.  This  section  uses  reduction  to  demonstrate  the 
existence  of  a  VPP  algorithm  without  actually  writing  out  the  algorithm.  As  a  consequence,  results  about 
C  are  generalized  to  C. 

Lemma  6  (C,AV{(jiJ^’~),u>^’~,Feas)  =vpp  (C jAV  Feas) 

Proof.  Let  11  =  {C,AV  ,u)^'=,Feas)  and  E'  =  {C',AV  jW^^^Feas). 

n'  PC  VPP  E:  C  C  C,  so  the  identity  mapping  is  a  reduction  from  E'  to  E. 

E  PC  VPP  E':  Let  I  =  (5, /,  wL=,  Af)  G  E  be  given,  where  S  is  defined  by  a  graph  G  =  {N,A,u,iJi,b)  G  C. 
Define  I'  =  (S' ,  f  ,  M)  G  E',  where  S'  is  defined  by  the  graph  G'  —  (N' ,  A' ,u' ,  p' ,b')  G  C.  The  graph 
G'  is  the  same  as  G,  except  as  follows: 

1.  (Eliminate  arcs  into  node  0)  For  each  arc  a  =  (v,  0)  G  A: 

(a)  Delete  arc  a  from  A' . 

(b)  Add  a  node  Wa  to  N'  with  b^^  =  0. 

(c)  Add  an  arc  a'  =  (i;,u;a)  to  A'  with  Ua'  =  Ua,  p.a<  =  A*a,  and  (xai)  =  /i,a 

2.  (Eliminate  parallel  arcs  from  node  0)  For  each  node  w  G  {1, . . . ,  m},  if  the  graph  contains  multiple  arcs 
from  node  0  to  node  v,  do  the  following  steps  for  all  but  one  of  the  arcs.  Let  a  =  (0,  v)  be  the  arc  on 
which  the  steps  are  being  performed. 

(a)  Delete  arc  a  from  A' . 

(b)  Add  a  node  Wa  to  N'  with  bw^  =  0. 

(c)  Add  an  arc  a'  =  (0,  Wa)  to  A'  with  «<,/  =  Ua,  Pa'  —  1.  and  (xa')  =  0. 

(d)  Add  an  arc  a"  =  {wa,  v)  to  A'  with  Ua"  —  Ua,  /ia"  =  Ma,  and  /(  ,,,,  (a?a")  =  fi,a  {xa")- 

3.  (Guarantee  an  arc  from  node  0  to  each  other  node)  For  each  node  v  £  N  such  that  arc  (0,  v)  ^  A:  Add 

an  arc  a'  =  (0,  u)  to  A!  with  Ua'  =  0,  fj,a'  =  L  and  /(  (xa’)  =  0. 

4.  (Guarantee  pairs  of  anti-parallel  arcs)  For  each  arc  {v,w)  G  A  such  that  arc  {w,v)  ^  A:  Add  an  arc 
a'  =  (w,  v)  to  A'  with  Ua'  —  0,  pa'  =  1)  and  (xq')  =  0. 

5.  (Guarantee  connectedness)  If  the  subgraph  of  G  induced  by  nodes  {1, . . . ,  n}  is  not  connected: 

(a)  Add  a  node  w  to  N'  with  bw  =  0. 
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(b)  Add  pairs  of  arcs  a!  —  {v,w)  and  a"  =  {w,v)  to  A',  with  Ua'  =  Wa”  =  0,  Ha'  =  Ha"  =  1,  and 
fi,a'  i^a')  =  /i,a»  (*o")  =  0,  to  at  least  one  node  v  in  each  component  of  the  subgraph. 

For  each  arc  a  that  is  in  both  A  and  A',  define  (xa)  =  /i,a  (®a);  then  /'  G  AV  By  construction, 

G'  €  C';  furthermore,  G'  can  be  constructed  in  VPP  time,  and  L(I)  is  polynomial  in  L(I').  The  reduction 
consists  of  constructing  G'  and  solving  T  as  described  in  Section  4.1.  This  yields  a  solution  of  I  directly 
using  the  correspondences  described  in  the  construction  of  G'.  ■ 

Theorem  3  and  Lemma  6  together  with  Lemma  4  in  [S094]  yield  the  following  corollary. 

CorollEiry  2  The  feasibility  problem  {C,A'P  ,  Teas)  has  a  VPP  algorithm  and  the  optimization 

problem  {C,AV  Opt)  has  an  FAS. 

Corollary  2  is  a  powerful  building  block.  An  explicit  VPP  algorithm  need  not  be  written  down  and  verified 
in  order  to  prove  the  existence  of  an  FAS  for  another  problem;  instead,  problems  of  interest  can  be  VPP 
reduced  to  the  problem  (C,AV  Feas).  This  is  done  for  several  STT  problems  in  Section  5. 

Several  knapsack,  job  scheduling,  and  lot-sizing  problems  are  treated  this  way  in  Sections  6  through  9. 


5  Easy  Problems  and  Hard  Problems 

Many  problems  can  be  defined  on  STT  networks  besides  those  discussed  in  the  previous  section.  This  section 
explores  the  boundary  between  “easy”  and  “hard”  problems.  The  easy  problems  are  those  for  which  an  FAS 
exists;  the  hard  problems  are  those  for  which  no  VPP  algorithm  can  exist  unless  V  =  JVV.  Section  5.1 
examines  problems  with  multiple  criteria  and  Section  5.2  considers  problems  in  which  the  flows  must  be 
exact.  A  summary  of  the  results  is  tabulated  in  Figure  6  at  the  end  of  this  section. 

The  proofs  of  hardness  use  the  Subset  Sum  problem,  which  is  A’P-complete[GJ79]. 

Subset  Sum 

Instance:  (n,  V,  v).  The  number  of  items  n;  the  target  sum  V;  and  the  sizes  Vj,  j  =  1, . .  .,n,  are  positive 
integers  such  that  vj  <  V ,  j  =  1, . . .  ,n. 

Question:  Find  a  set  S'  C  {1, . . . ,  n}  such  that  ^  vj  —  V. 

ies 


5.1  Multi- Criteria  Objectives 

This  section  extends  the  previous  results  about  C  to  the  multi-criteria  case.  Theorem  4  generalizes  Corollary  2 
to  include  multiple  criteria  and  order-reversing  functions.  Corollary  3  uses  Theorem  5  to  achieve  similar 
results  for  problems  with  deficient  flows.  Theorem  6  illustrates  that  problems  with  parallel  arcs  can  be  hard. 
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Theorem  4  For  each  r  €  and  w  €  fl’’,  the  problem  (C,AM.  (w)  ,w,  Opt)  has  an  FAS. 

Proof.  This  follows  from  applying  Lemma  4  and  Lemma  2  in  [S094]  to  Corollary  2  and  Theorem  2.  ■ 

A  similar  theorem  is  true  for  problems  with  deficient  flows;  it  will  be  proved  as  a  corollary  to  the  following 
theorem,  which  describes  a  more  general  relationship  between  problems  with  sufficient  flows  and  those  with 
deficient  flows. 

Theorem  5  Let  A*  €  {A1,AP},  C*  G  {Cl,CG},  and  M*  G  {Ml,  MG}.  Define  STT  networks  ^  = 
(A* /C* /M* /NM/FS)  and  =  {A* /C* / M*  j N M / F D).  For  any  direction  vector  u  G  LF ,  define  the 
problems  11  =  (w)  ,w,  Feas)  and  11'  =  (^',v4A4  (w)  ,  w,  Feas).  Then  11  =vppn'. 

Proof.  The  case  of  11  ocypp  H'  will  be  proved;  the  converse  is  proved  similarly.  Suppose  that  an  instance 
I  =  (S,  f,0J,  M)  G  n  is  given,  where  S  is  defined  by  the  graph  G  =  (N,  A,u,  fi,b).  Define  the  graph 
G'  =  {N,  A,  u,  p,  &'),  where  for  v  G  A  —  {0}, 

by  —  ^  ^  by  , 

a€A:  agA: 

a=(w,v)  a=(v,w) 

and  let  S'  be  the  feasible  set  defined  by  G'. 

Define  the  function  /'  by  fi{x)  —  fi{u  —  x),  i  =  1, ...  ,r.  Since  /  G  AM.  (w),  Lemmas  2  and  1  can  be  used 
to  show  that  /'  G  AM  (w).  So  I'  =  (S',  f',u},M)  G  H'. 

Since  f  (x)  =  f{u  —  x),  the  proof  will  be  completed  by  showing  that  a  flow  x  is  feasible  for  I  if  and  only  if 
the  flow  (u  —  x)  is  feasible  for  I'.  A  flow  x  is  feEisible  for  7  if  and  only  if  0  <  a:  <  u  and 

by  ^  0  .  (5) 

a^A:  agA: 

a={w,v)  at=:[v,w) 

But  —  u  ~  X  satisfies  0  <  x'  <  and 

^  ^  f^a  (^a  ^a)  ^  ^  (^a  ^a)  ~  ^ 

a^A:  a^A: 

a=(t7;,v)  a=(v,u') 

if  and  only  if  Equation  5  is  true.  This  demonstrates  the  claim,  so  n  ocvpp  H'.  ■ 

Note  that  the  construction  of  b'  requires  that  both  supply  and  demand  nodes  be  allowed,  so  the  fourth 
component  of  the  subclass  definition  in  Theorem  5  is  NM.  The  following  corollary  is  the  analog  of  Theorem  4 
for  problems  with  deficient  flows.  It  follows  from  applying  Lemma  3  and  Lemma  4  in  [S094]  to  Theorem  5 
and  Theorem  4. 

Corollary  3  Let  Ci  =  (Al/CG/MG/NM/FD).  Then  for  each  r  G  2+  and  u  G  LF ,  the  optimization 
problem  {Ci,AM  (w) ,  w.  Opt)  has  an  FAS. 
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Theorem  4  and  Corollary  3  cannot,  however,  be  extended  to  cover  STT  networks  that  contain  parallel  arcs, 
as  shown  by  the  following  theorem. 

Theorem  6  Let  Ci  =  {AP/CG/MlINQI FS)  and  Ca  =  {AP/CG/M\/NQ/FD).  Neither  the  problem 
(CijAV  Feas)  nor  the  problem  {C2,AP  ,  Feas)  has  a  VPP  algorithm  unless  V  = 

MV. 

Proof.  The  theorem  will  be  proved  by  reducing  Subset  Sum  to  each  of  the  STT  problems  in  such  a  way 
that  a  VPP  algorithm  for  one  of  the  latter  would  yield  a  polynomial-time  algorithm  for  the  former.  Since 
Subset  Sum  in  .AfP-complete,  such  a  VPP  algorithm  cannot  exist  unless  V  —  MV. 

First  consider  the  problem  Hi  =  {Ci,AV  ,w^>~,Feas).  Given  an  instance  (n,  V,  u)  of  Subset  Sum, 

construct  an  STT  network  G  =  {N,  A,u,  p,,b)  as  follows.  The  graph  is  pictured  in  Figure  3.  The  set  N  has 
three  nodes,  and  both  non-source  nodes  are  transshipment  nodes;  that  is,  N  =  {0,1,2},  with  6^  =  62  =  0. 
The  arcs  are  described  in  the  following  table: 


j 

Arc 

Capacity 

Multiplier 

Value  fj  {xj) 

1, . .  .,n 

(1,2) 

Vj 

1 

2n  -b  2  if  Xj  —  Vj 

2  if  0  <  Xj  <  Vj 

0  if  Xj  =  0 

n  -b  1 

(0,1) 

V 

1 

0 

n  +  2 

(2,0) 

V 

1 

1  if  Xj  =  V 

0  o.w. 

The  objective  function  /  (x)  =  fj  (xj)  is  in  AV  It  will  be  shown  that  the  instance  of  Subset 

Sum  has  a  solution  if  and  only  if  at  least  one  instance  (G,  k{2n  -|-  2)  -|-  l),  =  0, . . . ,  n,  of  Hi  has 

a  solution,  and  that  the  solution  to  any  of  these  instances  immediately  yields  a  solution  to  the  instance  of 
Subset  Sum.  The  reduction  of  the  instance  of  Subset  Sum  to  the  instance  of  Hi  is  polynomial  in  the  size  of 
the  Subset  Sum  instance,  so  a  VPP  algorithm  for  Hi  yields  a  polynomial-time  algorithm  for  Subset  Sum. 

If  the  Subset  Sum  instance  has  a  solution:  Suppose  that  S  C  (1, . . . ,  n},  with  |5|  =  fe,  is  a  solution  to  the 
instance  of  Subset  Sum.  One  feasible  solution  to  the  constructed  instance  of  Hi  with  k  —  k  is  xq^i  =  V, 
a^2,o  =  V,  Xj  =  Vj  for  j  S  S,  and  xj  =  0  for  y  G  (1, . . . ,  n}  —  S.  This  solution  has  value  (k{2n  +  2)  -f  1). 

If  the  instance  of  Hi  has  a  solution:  Suppose  that  a;  is  a  solution  to  the  instance  of  Hi  with  value  (k{2n  + 
2)-f  1).  Since  this  value  is  odd,  the  flow  in  arc  (2, 0)  must  be  V ;  this  leaves  k{2n  +  2)  of  the  objective  function 
unaccounted  for. 

Since  the  flow  into  node  2  is  sufficient,  at  least  V  units  of  flow  must  enter  node  2  on  arcs  1  through  n.  Since 
the  flow  into  node  1  is  sufficient,  and  since  uo,i  =  V,  the  total  flow  on  arcs  1  through  n  must  be  exactly  V. 
Since  the  remaining  part  of  the  objective  function  is  divisible  by  (2n-b2),  each  of  the  arcs  1  through  n  with 
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non-zero  flow  must  be  used  to  capacity;  no  subset  of  partially-full  arcs  could  contribute  a  factor  of  (2nd-2) 
to  the  objective  function.  The  indices  of  the  arcs  1  through  n  with  non-zero  flow  constitute  a  solution  of 
size  k  to  the  instance  of  Subset  Sum. 

The  second  part  of  the  theorem  is  about  the  problem  (C2,AV  (w'’~)  ,w^’~,Feas).  The  proof  is  similar, 
except  that  the  objective  functions  on  arcs  n  -f  1  and  n  -f  2  are  interchanged.  ■ 


5.2  Exact  Flows 

The  graphs  discussed  so  far  must  have  either  sufficient  or  deficient  flows.  This  section  considers  problems  in 
which  the  flows  must  be  exact,  that  is,  the  excess  flow  into  each  node  must  be  zero.  Theorem  7  identifies 
some  such  problems  that  can  be  solved  by  an  FAS.  Theorems  8  and  9  show  that  Theorem  7  cannot  be 
extended  to  graphs  with  both  supply  and  demand  nodes  or  to  problems  with  objective  functions  that  are 
not  order-preserving. 

Theorem  7  Lei  Ci  =  {Al/CG/Ml/N D/ FE)  and  C2  =  (Al/CG/Ml/NS/FE).  For  any  direction  vector 
w,  an  FAS  exists  for  each  of  the  problems  {Ci,AV  (w)  ,w,  Opt)  and  {C2,A'P  (w)  ,to,  Opt). 

Proof.  First  consider  the  problem  Hi  =  (Ci ,  AP  (w) ,  w.  Opt).  Let  r  be  the  number  of  components  in  w.  De¬ 
fine  the  STT  network  feasible  set  C3  =  (Al/GG/Ml/ND/FS)  and  the  problem  Ila  =  (C3,  AV  (u>)  ,u>,  Feas). 
Note  that  C3  C  C,  so  IIs  has  an  FAS. 

The  reduction  IIi  oc yppHs  will  be  proved.  By  Lemma4  in  [S094],  this  will  mean  that  Bi  has  an  FAS,  thereby 
proving  the  first  half  of  the  theorem.  The  proof  for  (C2,AV  (w)  ,ui,  Opt)  is  substantially  the  same,  and  so 
will  not  be  presented;  it  uses  a  reduction  to  a  problem  defined  on  the  feasible  set  {Al/GG /M\/N S/ FD). 

To  prove  the  reduction  Hi  ocyppHs,  suppose  that  an  instance  Ii  =  (^i,  /,  w)  G  Hi  is  specified,  where  Si  G  Ci 
is  defined  by  a  graph  G.  Let  S3  G  C3  also  be  defined  by  G,  and  define  the  instance  I3  =  {S3,  f,u,  M)  G  IIs 
for  any  M  G  2'''^ . 

Consider  any  flow  x  that  solves  I3.  This  flow  can  be  decomposed  into  flows  around  circuits  and  flows  on 
chains  (simple  directed  paths)  from  supply  nodes  to  demand  nodes[FF62,  AM093].  Suppose  that  the  flow 
around  some  circuit  in  such  a  decomposition  is  non-zero;  since  /  is  order-preserving  w.r.t.  w,  no  component 
of  the  objective  would  be  made  worse  by  eliminating  the  flow  around  the  circuit.  So  the  circuits  in  the 
decomposition  can  be  assumed  to  have  zero  flow. 

All  nodes  except  node  0  are  demand  nodes,  so  all  chains  in  the  decomposition  are  from  node  0  to  another 
node.  For  each  node  v  with  positive  excess,  that  is,  (*)  >  0,  consider  the  chains  from  node  0  to  node  v. 
Reducing  the  flow  on  one  or  more  of  these  chains  yields  a  flow  x'  with  e^  (x')  =  0.  As  with  the  circuits,  no 
component  of  the  objective  function  will  be  made  worse  by  reducing  the  flow  in  this  way. 
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Reduce  the  flow  in  this  way  for  each  node  to  obtain  an  exact  flow  x.  The  flow  x  solves  Ji,  so  Hi  ocvpp  Ha.  ■ 

The  next  theorem  shows  that  requiring  an  exact  flow  in  a  network  that  contains  both  supply  and  demand 
nodes  results  in  a  hard  problem. 


Theorem  8  Let  C\  =  {Al/CG/Ml/NM/FE)  he  an  STT  network.  Then  no  VPP  algorithm  exists  for  the 
problem  [Ci,AV  Peas')  unless  V  =  NV ■ 


Proof.  The  theorem  will  be  proved  by  reducing  Subset  Sum  to  the  problem  11  =  (Ci ,  AV  ,  Feas) 

in  such  a  way  that  a  VPP  algorithm  for  11  would  directly  yield  a  polynomial-time  algorithm  for  Subset  Sum. 
Since  the  latter  problem  is  A^'P-complete,  however,  such  a  VPP  algorithm  cannot  exist  unless  V  =  AfV. 

Given  an  instance  (n,  V,  v)  of  Subset  Sum,  an  STT  network  is  constructed,  as  illustrated  in  Figure  4.  Intu¬ 
itively,  the  constructed  network  has  a  node  with  supply  vj  for  each  item  j.  An  (n  -|-  l)-st  node  with  demand 
V  represents  the  desired  sum.  The  arcs  used  to  satisfy  the  demand  define  a  solution  to  the  instance  of  Subset 
Sum. 

More  formally,  construct  an  STT  network  G  —  {N,A,u,fi,b)  as  follows.  The  set  N  has  n  +  2  nodes 
{0, 1, . . . ,  n,  n-h  1}.  Node  j  has  supply  vj,  that  is,  bj  —  —vj,  j  —  1, . . . ,  n.  Node  n  -f- 1  has  demand  V,  that 
is,  6„+i  =  V.  The  set  A  has  2n  arcs;  they  are  aj  =  (j,  0)  and  /3j  =  (i,  n  +  1)  for  j  =  1, . . . ,  n.  Arcs  aj  and 
fij  have  capacity  vj,  multiplier  1,  and  the  following  objective  function; 


9{xk)  = 


if  a;*,  >  1 
if  a:*  =  0  . 


The  objective  function  /  (a:)  =  [s'  +  9  is  in  AV  It  will  be  shown  that  the  instance 

of  Subset  Sum  has  a  solution  if  and  only  if  the  instance  (G, /,  n)  of  11  has  a  solution,  and  that  the 
solution  to  this  instance  immediately  yields  the  solution  to  the  instance  of  Subset  Sum.  The  reduction  of 
the  instance  of  Subset  Sum  to  the  instance  of  11  is  polynomial  in  the  size  of  the  Subset  Sum  instance,  so  a 
VPP  algorithm  for  11  yields  a  polynomial-time  algorithm  for  Subset  Sum. 

If  the  Subset  Sum  instance  has  a  solution:  Suppose  that  S  C  {1, . . . ,  n}  is  a  solution  to  the  instance  of  Subset 
Sum.  One  feasible  solution  to  the  constructed  instance  of  11  is  xp.  =  vj  for  j  G  S,  xp.  =  0  for  j  ^  S,  and 
Xaj  =  Vj  —  xp- ,  j  =  1, . . . ,  n.  This  solution  has  value  n. 

If  the  instance  of  11  has  a  solution:  Suppose  that  a;  is  a  solution  to  the  instance  of  11  with  value  not  greater 
than  n.  Each  node  j  must  dispose  of  its  supply  on  arc  aj ,  ,  or  both.  The  objective  function  must  therefore 

include  at  least  unit  value  for  the  flow  leaving  each  of  these  nodes,  and  hence  must  be  at  least  n. 

The  solution  to  11  therefore  has  value  exactly  n,  which  means  that  each  of  the  nodes  1, . .  .,n  uses  exactly 
one  of  its  outgoing  arcs.  Let  S  be  the  index  set  of  the  arcs  /3j  that  have  non-zero  flow.  Then  for  j  G  S,  it 
must  be  the  case  that  xpj  =  Vj,  and  so  =  V .  So  S  constitutes  a  solution  to  the  instance  of  Subset 
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Figure  3:  The  graph  constructed  in  the  proof  of  Theorem  6. 


Figure  4:  The  graph  constructed  in  the  proof  of  Theorem  8. 
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Sum.  ■ 

Similarly,  the  next  theorem  shows  that  using  objective  functions  that  are  not  order-preserving  on  networks 
that  contain  only  transshipment  nodes  results  in  a  hard  problem. 

Theorem  9  Let  Ci  =  (Al/CG/Ml/NQ/FE)  be  an  STT  network.  Then  no  VPP  algorithm  exists  for  the 
problem  (Ci^AV  unless  V  =MV. 

Proof.  The  theorem  will  be  proved  by  reducing  Subset  Sum  to  the  problem  11  =  (Ci ,  AV  (wb<)  ,wb>  ,  Feas) 
in  such  a  way  that  a  VPP  algorithm  for  11  would  directly  yield  a  polynomial-time  algorithm  for  Subset  Sum. 
Since  the  latter  problem  is  A/’P-complete,  however,  such  a  VPP  algorithm  cannot  exist  unless  V  —  MV. 

Given  an  instance  (n,  V,  v)  of  Subset  Sum,  an  STT  network  is  constructed,  as  pictured  in  Figure  5.  Intuitively, 
each  node  j  receives  flow  Vj  from  node  0  and  either  sends  it  on  to  node  (n  +  1)  or  returns  it  to  node  0.  The 
arcs  into  node  n  H- 1  correspond  to  a  solution  for  the  instance  of  Subset  Sum. 

More  formally,  construct  an  STT  network  G  =  {N,A,u,n,b)  as  follows.  The  set  N  has  n  -f  2  nodes 
{0, 1, . . . ,  n,  n  -f  1}.  All  nodes  are  transshipment  nodes,  that  is,  bj  =  0,  j  =  1, ...  ,n  +  1.  The  set  A  has 
3n  -b  1  arcs  that  are  described  in  the  following  table: 


Arc 

Capacity 

Multiplier 

Value  fk  (xk) 

(0,  j):  i  =  l,...,n 

1 

0 

(j,  0),  j  =  l,...,n 

U,n  +  1),  j  =  l,...,n 

1 

1  if  Xik  =  Vj 

0  o.w. 

(n  -b  1,  0) 

V 

1 

1  iixk  =  V 

0  o.w. 

The  objective  function  /  (a;)  =  fk  (®fc)  in  •  I*  will  be  shown  that  the  instance  of  Subset 

Sum  has  a  solution  if  and  only  if  the  instance  (G,  n  +  l)  of  H  has  a  solution,  and  that  the  solution  to 

this  instance  immediately  yields  the  solution  to  the  instance  of  Subset  Sum.  The  reduction  of  the  instance 
of  Subset  Sum  to  the  instance  of  11  is  polynomial  in  the  size  of  the  Subset  Sum  instance,  so  a  VPP  algorithm 
for  n  yields  a  polynomial-time  algorithm  for  Subset  Sum. 

If  the  Subset  Sum  instance  has  a  solution:  Suppose  that  SC  {l,...,n}isa  solution  to  the  instance  of  Subset 
Sum.  One  feasible  solution  to  the  constructed  instance  of  11  is  Kj.n+i  =  Vj  for  j  E  S  and  Xj^n+i  —  0  for 
j  0  S;  —  V;  and  for  j  =  1, . .  .,n,  xqj  =  Vj  and  xj^o  =  vj  —  xj^n+i-  This  solution  has  value  n  -|-  1, 

and  so  solves  the  instance  of  11. 

If  the  instance  of  11  has  a  solution:  Suppose  that  a;  is  a  solution  to  the  instance  of  n  with  value  at  least  n  +  l. 
It  must  be  the  case  that  o  =  V  and  that  each  node  j,  j  =  I,  .. .,  n,  sends  flow  vj  in  exactly  one  of  the 
arcs  (j,  n  -b  1)  and  {j,  0)  and  zero  flow  in  the  other.  In  order  for  this  to  happen,  the  flow  in  each  arc  (0,  j), 
j  =  1,  . . .,  n,  must  be  vj . 
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Figure  5:  The  graph  constructed  in  the  proof  of  Theorem  9. 
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Cl 

n 

Difficulty 

Reference 

( A1  /CG/MG/NM/FS) 
(Al/GG/MG/NM/FD) 
(AP/GG/Ml/m/FS) 
(AP/GG/Ml/NO/FD) 

(Cl ,  AM  (w’’’-) ,  w'’-- ,  Opt) 
(Ci,^A4  (a;),w,Opt) 
(Ci,.4P(a;L=i  a;L=,Feas' 
{Ci,AV  Feas 

Easy 

Easy 

Hard 

Hard 

Theorem  4 
Corollary  3 
Theorem  6 
Theorem  6 

(AlIGGfMllNDlFE) 

(Al/GG/Ml/NS/FE) 

(Al/GG/Ml/NM/FE) 

(Al/GG/Ml/NQ/FE) 

{Ci,AP (w)  ,w,Opt) 
{Ci,AV{oj)  ,w,Opt) 
(Ci,>lP(wb<)_a;b<,Feas' 
(Ci.^T’  (wb<)  ,wi.>_Feas^ 

Easy 

Easy 

Hard 

Hard 

Theorem  7 
Theorem  7 
Theorem  8 
Theorem  9 

Figure  6:  Summary  of  STT  complexity  results. 


Let  S  be  the  set  of  nodes  j,  j  =  1, . . . ,  n,  that  send  non-zero  flow  in  arc  {j,  n  -)-  1).  Since  6„+i  =  0  and  the 
excess  at  each  node  must  be  0,  it  follows  that  so  S'  constitutes  a  solution  to  the  instance  of 

Subset  Sum.  ■ 

The  results  of  this  section  are  summarized  in  Figure  6. 


6  Knapsack  Covering  Problems 

This  section  discusses  the  form  of  constraints  that  deflne  the  feasible  regions  for  knapsack  problems,  as  a 
prelude  to  the  discussion  of  specific  problems  in  the  following  two  sections.  The  knapsack  problems  considered 
in  this  paper  are  special  cases  of  the  Multi-Dimensional  Knapsack  problem,  which  is  formulated  below.  In 
general,  no  FAS  can  exist  for  this  problem  unless  P  =  NV . 


Integer  Multi-Dimensional  Knapsack  Feasible  Set 

Instance:  (n,  m,  V,  v,  u).  The  number  of  items  n  and  number  of  constraints  m  are  non-negative  integers. 

The  constraint  upper  bounds  Vi,  i  =  1, . . . ,  m;  the  item  values  Vij ,  i  =  1, . . .  ,m.,  j  —  1, . . .  ,n', 
and  the  item  quantity  upper  bounds  uj,  j  =  l,...,n;  are  non-negative  integers  such  that 
Vij  <Vi,i=l,...,m,j=l,...,n. 

Feasible  Set:  S{n,  m,  V,  v,  u)  is  the  set  of  all  x  E  2^"“^  that  satisfy  the  following  constraints: 


n 


E 


< 

. ,  m 

< 

Uj, 

II 

..,n 

This  generic  knapsack  problem  has  “<”  inequality  constraints.  This  is  characteristic  of  the  “packing”  form  of 
knapsack  problems,  in  that  upper  limits  are  placed  on  the  amounts  that  can  be  packed  into  a  knapsack.  Many 
interesting  problems,  however,  have  feasible  regions  that  are  naturally  expressed  by  knapsack  constraints 
with  “>”  inequalities;  these  are  called  knapsack  “covering”  problems. 
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Definition  7  Let  ^  be  a  family  of  feasible  regions  such  that  each  5  €  ^  can  be  described  as  the  set  of  all 
X  €  2’”+  that  satisfy,  for  some  integers  n,  mi,  m2;  Vi,  i  =  l,...,mi;  Vij,  i  =  l,...,mi,  j  =  l,...,n;  VFi, 
i  =  1, . . m2;  Wij,  i  —  I,..  .,m2,  j  =  1, . . . ,  n;  and  Uj ,  j  =  1, . . . ,  n; 

n 

VijXj  <  V{,  i  =  1, . . .  ,mi 

j=i 

WijXj  =  Wi,  i=l,...,m2 
t=i 

Xj  <  Uj,  j  =l,...,n  . 

Such  a  set  S  is  said  to  be  of  the  packing  type.  The  covering  form  of  S,  written  is  the  feasible  region 
defined  as  the  set  of  all  j/  S  that  satisfy 


n  n 


E 

Vijyj 

> 

Y^VijUj 

-Vi,  i  = 

.,mi 

i=i 

n 

E 

WijVj 

— 

i=i 

n 

Y^WijUj 

1 

II 

.,m2 

j  =  l 

Vi 

< 

3  =  1 

.  .,n  , 

and  the  covering  form  of  ^  is 

^cov  _  I  ^cov  .  5  g  ^  j. 


The  present  work  primarily  considers  knapsack  problems  of  the  packing  form.  The  following  theorem  shows 
that  this  focus  typically  does  not  sacrifice  much  generality. 

Theorem  10  Let  ^  be  a  family  of  feasible  regions  of  the  packing  type  and  w  G  0-  Then 

(1)  i^,AV{co),ij,Opt)  =VPP 

(2)  i^,ATl{w),u,Opf)  =VPP  {^^°\AV{w),u;,Opf). 

(3)  (^,AMiu>),u,Opt)  =VPP  {^^°\AM{w),Lo,Opt). 

Proof.  Statement  (3)  follows  directly  from  statements  (1)  and  (2).  The  proof  that  f^,AV  (w)  ,U),  Opt)  oc 
vpp  ,  ATZ  (u>)  ,u,  Opt)  is  given  here;  the  remaining  three  reductions  in  (1)  and  (2)  are  proved  similarly. 

Let  (S,  /,  w)  G  ,  AV  (w) ,  w,  Opt).  For  any  a:  G  5,  let  j/  =  u  —  a;.  Then  y  G  ;  for  example, 

n  n 

=  Yl'VijiUj  -  Xj) 

i=i  j=i 

n 

i=i 
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Now  define  the  function  f(y)  —  f  [u  —  y).  Then  f  E  ATZ  (w):  if  j/i ,  2/2  €  5*^°^ ,  then 

2/1  <  2/2  =>  u-y2<u-yi 

Du  (/  {u  -  2/2) ,  /  (u  -  2/1))  because  /  6  AT  (w) 

Duif'{y2),f'{yi))  • 

Furthermore,  f{u  —  x)  =  f{x).  So  (^, (w) , w, Opt)  aypp  ,w,Opt).  ■ 

Many  kinds  of  problems  admit  an  FAS  if  the  feasible  region  has  only  packing  constraints  or  only  covering 
constraints.  Combining  the  two  kinds  of  constraints,  however,  often  yields  problems  that  have  no  FAS  unless 
T  =  AfV.  That  is  why  knapsack  problems  in  which  the  knapsack  constraint  is  an  equality  are  generally 
difficult  to  solve. 

For  some  problems,  such  as  the  Partially  Ordered  Knapsack  problem  discussed  in  Section  8,  not  all  the 
constraints  are  reversed  to  define  the  covering  form  of  the  problem.  Theorem  10  cannot  be  directly  applied 
to  such  problems,  but  it  can  sometimes  be  used  indirectly,  as  it  is  in  Section  8. 


7  The  Arborescent  Knapsack  Problem 

This  section  and  the  next  discuss  the  existence  and  non-existence  of  FAS’s  for  a  variety  of  knapsack  and 
scheduling  problems.  Note  that  the  results  presented  are  existential.  The  characteristics  of  a  particular 
problem  must  be  considered  carefully  in  order  to  develop  a  practical  algorithm  for  its  solution.  An  excellent 
example  of  such  an  analysis  appears  in  [Law79]. 

In  the  Arborescent  Knapsack  problem,  items  are  grouped  into  sub-knapsacks,  or  stutfsacks.  Stuffsacks  may 
be  nested  within  other  stuffsacks  and  each  stuffsack  has  its  own  capacity  bound.  This  structure  is  called 
“arborescent”  because  the  “nested  within”  relationship  on  the  stuffsacks  can  be  represented  by  a  directed  tree, 
or  arborescence,  as  illustrated  in  Example  2  later  in  this  section.  After  defining  the  Arborescent  Knapsack 
problem  and  proving  that  it  has  an  FAS,  two  special  cases  that  arise  in  the  context  of  job  scheduling  are 
discussed. 


7.1  The  General  Problem 

After  formalizing  the  notion  of  stuffsacks,  the  Integer  Arborescent  Knapsack  feasible  region  is  defined  and 
the  existence  of  FAS’s  for  certain  problems  defined  on  it  is  proved. 

Definition  8  A  collection  of  sets  B  =  {Bi, . . . ,  Bm}  is  called  arborescent  if  for  each  pair  i,j  E  {1, . . . ,  m}, 
with  i  ^  j,  one  of  the  following  relations  holds:  Bi  C  Bj,  Bj  C  Bi,  or  Bi  D  5^  =  0.  B  is  sometimes  called  a 
collection  of  stuffsacks.  ■ 

The  Arborescent  Knapsack  feasible  set  may  now  be  defined  formally  as  follows: 
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Integer  Arborescent  Knapsack  Feasible  Set 

Instance:  (n,  m,  B,  V,  v,  u).  The  number  of  items  n  and  number  of  stuffsacks  m  are  non-negative  integers. 

The  collection  of  stuffsacks  B  is  an  arborescent  collection  of  m  sets  5,-  C  n},  i  = 

1, . . . ,  m,  with  Bm  =  n}.  The  stuffsack  volumes  Vi,  i  —  I, . . m;  the  item  volumes  Vj, 

j  =  1, . . . ,  n;  and  the  item  quantity  upper  bounds  uj,  j  =  1, . , . ,  n;  are  non-negative  integers 
such  that  if  j  €  Bi,  then  Vj  <  Vi,  «  =  1,  . . m,  j  —  1,  . . n. 

Feasible  Set:  S{n,  m,  B,  V,  v,  u)  is  the  set  of  all  x  £  that  satisfy  the  following  constraints: 

VjXj  <  Vi,  i  =  l,...,m 

j€Bi 

<  Uj,  j  =  l,...,n. 

The  nesting  property  of  the  sets  Bj  gives  the  multiple  constraints  a  special  structure  that  allows  the  Arbores¬ 
cent  Knapsack  problem  to  have  an  FAS,  despite  the  difficulty  of  the  Multi-Dimensional  Knapsack  problem. 
An  FAS  that  uses  0{n^/e^)  time  for  the  binary  version  of  this  problem  with  a  single  linear,  additively 
separable  criterion  is  presented  in  [GL79b].  The  following  theorem  generalizes  their  result  to  the  integer 
case,  in  which  an  item  may  appear  more  than  once  in  the  knapsack  (up  to  Uj  copies  of  item  j  may  he  used). 
In  addition,  this  theorem  accommodates  multiple  criteria. 

Theorem  11  Let  ^  be  the  collection  of  Integer  Arborescent  Knapsack  feasible  sets  and  w  E  fi.  Then 
{'^,AA4  (w)  ,u>,  Opt)  has  an  FAS. 

Proof.  Let  Hi  =  (w) ,  w,  Feas)  and  II2  =  (C2,  MAI  (w)  ,  w,  Feas),  where  the  feasible  set  C2  = 

(Al/CG/MG/NO/FD).  The  reduction  Hi  ocvppn2  will  be  shown.  By  Lemma4  in  [S094]  and  Corollary  3, 
this  reduction  will  prove  the  present  theorem. 

Suppose  that  an  instance  B  -  {S,f,u),M)  of  IIi  is  given,  with  S  =  S {n,m,  B,V,v,u).  Assume  that  the 
sets  Bi  are  “topologically”  ordered,  so  that  B,-  C  Bk  implies  that  i  <  k]  such  an  ordering  can  be  found 
quickly[AM093].  Define  the  “first  containing  set”  function  <j)  for  items  by 

(t>{j)  =  min{  :  j  €  Bk)  ,  j  -  1, . .  .,n 

and  a  similar  function  cr  for  sets  by 

I  min  {k  :  Bi  C  Bk}  i  =  I, . . . ,  m  -  1 

a{t)  =  I 

0  t  =:  m  . 


Define  an  instance  I2  =  {G2,  f  ,oj,M)  of  112,  where  G2  is  defined  by  an  STT  network  G  =  {N ,  A,u,  p.,b). 
An  example  of  the  following  construction  is  given  in  Example  2.  The  set  N  has  (n  -|-  m  4-1)  nodes:  a  node 
aj  corresponding  to  each  item  j,  a  node  /?,  corresponding  to  each  stuffsack  i,  and  node  0.  All  nodes  are 
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transshipment  nodes.  The  set  A  has  (2n  +  m)  arcs  that  are  described  in  the  following  table: 


Arc 

Capacity 

Multiplier 

Value  (a;*,) 

(0,a^), 

i  =  1,  • 

.,n 

Uj 

fs  i^j) 

{Oj  , 

j  =  1,  • 

.,n 

VjUj 

1 

0 

(A’j 

i  =  1, . . 

, ,  m 

Vi 

1 

0 

If  h  has  a  solution:  Suppose  that  x  is  a  feasible  solution  to  Ii .  One  feasible  solution  to  I2  is  shown  in  the 
following  table:  _ _ 


Arc 

Flow 

(0,  cvj)j 

i  =  1, •  ■ 

.,n 

Xj 

[aj , 

j  • 

.,n 

VjXj 

(A’ ) 

?■=!,.. 

. ,  m 

Y 

j€Bi 

This  is  a  feasible  flow  that  has  the  same  objective  function  value  as  does  the  solution  of  Ii. 


If  I2  has  a  solution:  Suppose  that  y  is  a  feasible  flow  for  I2 ,  and  consider  the  following  assignments  for  the 
knapsack  variables: 

j  =  1,  .  .  ■  (7) 

The  capacity  constraints  on  arc  (0, aj)  guarantee  that  Xj  <  Uj,  j  =  .,n.  It  will  be  shown,  by  induction, 

that  Ylj^B  ‘^i^j  —  ypiA<T(,i)>  *  ~  1, . . . ,  m.  Then,  by  the  capacity  constraints  on  the  arcs  (/?,-, it  follows 
that  VjXj  <  Vi,  so  a;  is  a  feasible  solution  for  h  with  the  same  objective  function  as  y  has  in  l2- 


i  =  1  :  VjXj 

j€Bi 


jeBi 

jSBi 

j6Bi 


by  Equation  7 
nodes  have  deficient  excess 
sets  are  topologically  ordered 
node  /?i  has  deficient  excess. 


i>l:  Y 

j€Bi 


— 

Y 

+ 

E 

Y 

arborescent  network  structure 

j:<f>(j)-i 

k:<7{k)zzi 

j€Bk 

< 

+ 

E 

Y 

nodes  have  deficient  excess 

3-4U)=i 

k:(7{k)=i 

j&Bh 

< 

Y  y^Bpi 

+ 

E 

Y 

by  induction 

j:<t>{j)=i 

k:<7{k)=i 

j^Bh 

< 

ypiA.M 

node  A  has  deficient  excess. 
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Example  2  (Network  construction  for  the  proof  of  Theorem  11)  Let  n  =  7  and  m  —  5,  with  the 
following  arborescent  collection  of  sets: 


Bi  =  {3,4,5} 

B2  =  {6} 

Bs  =  {1,2} 

54  =  {1,2, 3, 4, 5} 

Sb  =  {1,2, 3, 4, 5, 6, 7}. 


Then  the  functions  and  cr(-)  are  defined  as  follows; 


Item  j 

Set 

Set  i 

Set  cT{i) 

1 

3 

1 

4 

2 

3 

2 

5 

3 

1 

3 

4 

4 

1 

4 

5 

5 

1 

5 

0 

6 

2 

7 

5 

and  the  corresponding  network  is  shown  in  Figure  7.  ■ 


7.2  The  Maximum  Job  Sequencing  with  Due  Dates  Problem 

A  special  case  of  the  Arborescent  Knapsack  problem  is  the  Triangular  Knapsack  problem,  in  which  the 
stuffsacks  are  nested  sequentially.  More  precisely,  the  Triangular  Knapsack  problem  is  an  Arborescent 
Knapsack  problem  in  which  m  —  n  and  Bi  —  {1, . .  .,i}.  One  reason  for  interest  in  this  problem  is  that  a 
common  scheduling  problem  can  be  formulated  as  a  Triangular  Knapsack  problem.  Details  about  related 
scheduling  problems  can  be  found  in  [GLLRK79]. 

In  the  Max  Job  Sequencing  with  Due  Dates  problem,  a  set  of  independent  jobs  must  be  scheduled  on  a 
single  machine.  Those  completed  by  their  due  dates  earn  a  profit,  but  those  completed  late  earn  nothing; 
the  objective  is  to  earn  as  much  profit  as  possible.  This  can  be  described  more  formally  as  follows: 

Max  Job  Sequencing  with  Due  Dates  (Max  JSD),  Version  1 

Instance:  (n,p,  t,  d).  The  number  of  jobs  n;  the  profits  pj,  j  —  1, . . . ,  n;  the  execution  times  tj,  j  =  1, . . . ,  n; 

and  the  due  dates  dj,  j  =  are  non-negative  integers  such  that  tj  <  dj ,  j  —  1, . . .  ,n,  and 

dj  ^  dj^i,  j  —  1 , . . . ,  (n  1). 

Question:  Schedule  the  jobs  to  earn  the  maximum  possible  profit,  where  job  j  earns  profit  pj  if  the  job  is 
completed  before  time  dj  and  earns  0  profit  otherwise. 
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Figure  7:  Network  constructed  as  in  the  proof  of  Theorem  11  for  the  instance  of  Example  2. 
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An  FAS  that  runs  in  time  O  (n^/e)  for  this  scheduling  problem  is  described  in  [Sah76,  Sah77]. 

To  see  that  this  problem  can  be  formulated  as  a  Triangular  Knapsack  problem,  recall  that  the  profit  is 
independent  of  the  order  in  which  the  jobs  are  performed,  and  depends  only  on  the  relationship  of  the  finish 
times  to  the  due  dates.  The  set  of  optimal  schedules  must  therefore  include  at  leeist  one  schedule  in  which 
all  the  jobs  that  are  completed  on  time  are  performed  before  those  that  are  finished  late;  so  the  search  for 
an  optimal  schedule  can  be  restricted  to  schedules  of  this  kind.  Define  the  variables  xj ,  j  =  1, . . .  ,n,  hy 

{1  if  job  j  is  completed  on  time, 

(8) 

0  otherwise. 

The  Max  JSD  problem  can  then  be  formulated  as: 

Max  Job  Sequencing  with  Due  Dates  (Max  JSD),  Version  2 

Instance:  {n,p,t,d).  The  number  of  jobs  n;  the  profits  p^,  j  =  1, . .  .,n;  the  execution  times  j  =  1, . .  .,n; 

and  the  due  dates  dj,  j  —  1, . . . ,  n;  are  non-negative  integers  such  that  tj  <  dj ,  j  =  1, . . .  ,n,  and 
dj  <  dj+i,  i  =  l,...,(n-  1). 

Question:  Find  x  to  solve  « 

max  E  PjXj 

j=l 

3 

S.t.  ^  dj  J  j  —  1,  ...  J  72 

ifc  =  l 

€  {0,1},  i=l,...,72 

This  has  the  form  of  a  Triangular  Knapsack  problem.  This  formulation  is  similar  to  that  alluded  to  in 
[Sah76,  Sah77,  GL79a]  and  described  more  fully  in  [Rot64,  LM69,  GL78,  GL81]  for  a  similar  scheduling 
problem  that  is  discussed  in  the  next  section.  The  results  of  Section  7.1  imply  the  existence  of  an  FAS  for 
the  multi-criteria  form  of  this  problem. 

7.3  The  Minimum  Job  Sequencing  with  Due  Dates  Problem 

The  Job  Sequencing  with  Due  Dates  Problem  can  be  approached  differently  by  assessing  penalties  for  jobs 
that  are  completed  late,  rather  than  providing  profits  for  jobs  that  are  completed  on  time.  The  objective  in 
this  case  is  to  minimize  the  sum  of  the  penalties,  and  the  problem  is  called  the  Min  Job  Sequencing  with 
Due  Dates  (Min  JSD)  problem.  An  exponential-time  dynamic  programming  algorithm  for  this  problem  is 
presented  in  [LM69].  An  FAS  that  requires  O  {n^/e)  time  is  described  in  [GL78];  better  algorithms  that  run 
in  time  O  (72^1og(n)  -|-  rfi  jt)  are  given  in  [GL79a,  GL81]. 

The  Min  JSD  problem  is  not  just  a  trivial  variant  of  the  Max  JSD  problem,  because  the  complexities  of 
approximating  the  solution  of  a  minimization  problem  and  of  the  corresponding  maximization  problem  can 
be  quite  different.  It  has  been  speculated  that  approximating  solutions  for  the  Min  JSD  problem  may  be 
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inherently  harder  than  doing  so  for  the  Max  JSD  problem[GL81].  Nevertheless,  the  approximation  schemes 
of  [GL78,  GL79a,  GL81]  are  derived  by  reducing  the  Min  JSD  problem  to  the  Min  Arborescent  Knapsack 
problem.  The  variables  are  the  complements  of  those  used  in  Equation  8: 


yj  = 


1  if  job  j  is  completed  late, 
0  otherwise. 


(9) 


The  problem  can  be  formulated  as  follows: 


Min  Job  Sequencing  with  Due  Dates  (Min  JSD),  Version  1 

Instance:  {n,p,t,d).  The  number  of  jobs  n;  the  penalties  pj,  j  =  l,...,n;  the  execution  times  tj,  j  = 
1, . . . ,  n;  and  the  due  dates  dj,  j  =  1, . . . ,  n;  are  non-negative  integers  such  that  tj  <  dj ,  j  = 
1, . . . ,  n,  and  dj  <  dj+i,  j  =  1, . . .,  (n  —  1). 

Question:  Find  y  to  solve  " 

min  E  PjVj 

3 

s.t.  ^  <  dj,  j  =  l,...,n 

k=l 

Vj  G  {0,1},  i  =  l,...,Tt, 

which  is  a  Triangular  Knapsack  problem.  Setting  dj  =  I2i=i  —  dj,  however,  allows  the  Min  JSD  problem 
to  be  formulated  as  a  knapsack  covering  problem: 


Min  Job  Sequencing  with  Due  Dates  (Min  JSD),  Version  2 

Instance:  {n,p,t,  d).  The  number  of  jobs  n;  the  profits  pj,  j  —  1, . . . ,  n;  the  execution  times  tj,  j  =  1, . . . ,  n; 

and  the  due  dates  dj,  j  =  1, . . . ,  n;  are  non-negative  integers  such  that  tj  <  dj,  j  =  1, . . . ,  n,  and 
dj  <dj+i,j-  l,...,(n-l). 

Question:  Find  y  to  solve  ^ 

min  E  PjVj 

3  =  1 
3 

S.t.  ^  ^j>  j  =  l,...,n 

k=l 

Vj  G  {0,1},  i  =  l,...,n. 

This  is  a  Binary  Arborescent  Knapsack  Covering  problem.  To  our  knowledge  the  Arborescent  Knapsack 
Covering  problem  has  only  been  treated  by  considering  the  special  case  of  this  scheduling  problem.  It  does, 
however,  satisfy  the  preconditions  of  Theorem  10.  The  previously  known  results  for  the  Min  JSD  problem 
can  therefore  be  generalized  to  guarantee  an  FAS  for  the  Integer  Arborescent  Knapsack  Covering  problem 
with  multiple  criteria. 
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8  The  Partially  Ordered  Knapsack  Problem 

Another  special  case  of  the  Multi-Dimensional  Knapsack  problem  whose  structure  allows  for  the  existence 
of  an  FAS  is  the  Partially  Ordered  Knapsack  problem.  This  is  a  variant  in  which  items  must  be  chosen  in 
accordance  with  specified  precedence  constraints.  Such  a  constraint  requires  that,  for  instance,  if  an  item  j 
is  selected,  then  another  specified  item  i  must  also  be  selected. 

Requirements  of  this  nature  are  quite  common,  and  can  arise  in  many  ways.  One  source  of  precedence 
constraints  is  physical  structure.  For  example,  a  manufacturer  must  decide  which  items  to  produce  in  a 
fixed  time  period,  but  must  produce  the  components  for  certain  complex  assemblies  before  producing  tbe 
assemblies  themselves. 

Precedence  constraints  define  a  partial  order  on  the  set  of  items;  if  item  i  must  be  selected  whenever 
item  j  is  selected,  then  i  j  (read  i  precedes  j).  Such  a  partial  order  can  be  represented  by  a  directed 
graph  G(~>)  =  (TV,  A);  TV  contains  a  node  for  each  item,  and  for  each  i,j  E  N,  A  contains  arc  {i,j)  if  and 
only  if  i  j.  The  graph  G{^)  is  assumed  to  be  acyclic,  that  is,  there  is  no  group  of  items,  all  of  which 
must  either  be  chosen  or  not  chosen.  The  graph  is  also  assumed  to  be  connected.  If  this  is  not  the  case, 
consider  the  maximal  connected  components  of  G('^);  a  “dummy”  item  can  be  added  that  precedes,  or  is 
preceded  by,  one  item  from  each  component. 

Several  kinds  of  partial  orders  are  particularly  interesting. 

Definition  9  A  partial  order  is  a  tree  partial  order  if  the  undirected  version  of  G('^)  is  a  tree.  ■ 

Definition  10  A  partial  order  is  an  in-tree  partial  order  if  it  is  a  tree  partial  order  and  G('^)  contains 
some  node  (the  “root”)  toward  which  all  its  arcs  point.  ■ 

Definition  11  A  partial  order  is  an  out-tree  partial  order  if  it  is  a  tree  partial  order  and  G(-^)  contains 
some  node  (the  “root”)  away  from  which  all  its  arcs  point.  ■ 

Figure  8  shows  an  in-tree  partial  order.  Figure  9  shows  an  out-tree  partial  order,  and  Figure  10  shows  a 
general  tree  partial  order  that  is  neither  an  in- tree  partial  order  nor  an  out-tree  partial  order.  Figure  11 
shows  a  partial  order  that  is  not  a  tree  partial  order. 

The  feasible  set  for  a  Partially  Ordered  Knapsack  problem  is  defined  as  follows. 

Partially  Ordered  Knapsack  Feasible  Set 

Instance:  5(n,IV,  v,^).  The  number  of  items  n;  the  knapsack  volume  P;  and  the  item  volumes  vj, 

j  =  1, . . . ,  n;  are  non-negative  integers  such  that  vj  <  V ,  j  =  1, . . . ,  n.  The  partial  order  ^  is 
defined  on  the  set  of  items  {1, . . . ,  n}. 
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Figure  8:  An  in- tree  partial  order. 


Figure  9:  An  out-tree  partial  order. 
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Feasible  Set:  S{n,  V,v,'^)  is  the  set  of  all  x  G  satisfying 

n 

j-l 

Xi  >  Xj  Zj  —  1, . . . ,  n  such  that  i  ^  j 

Xj  G  {0,1}  i=l,...,n 

The  Partially  Ordered  Knapsack  problem  with  a  single  criterion  of  the  form  J2j-iPjXj  is  strongly  TV'P-hard 
when  general  partial  orders  are  allowed,  but  can  be  solved  in  pseudo-polynomial  time  when  only  tree  partial 
orders  are  allowed [GJ 79].  A  PTAS  for  the  maximization  form  of  this  single-criterion  problem,  restricted  to 
in-tree  partial  orders,  is  described  in  [IK78];  it  runs  in  time  O  The  algorithm  of  [IK78]  is  also 

applied  to  the  minimization  form  of  the  problem  in  which  the  knapsack  constraint  is  of  the  covering  type 
and  the  precedence  constraints  form  an  in-tree  partial  order. 

A  variety  of  stronger  results  is  presented  in  [JN83].  An  FAS  is  presented  for  out-tree  or  in-tree  partial 
orders.  The  FAS  is  based  on  the  standard  “bottom-up”  dynamic  programming  approach  and  runs  in  time 
0{n^/€^  +  n^log{P*)).  A  new  dynamic  programming  paradigm,  the  “left-right”  approach,  that  yields  an 
FAS  that  runs  in  time  O  (n^/e^ -f  log(P*))  for  out-tree  partial  orders  and  0(n^/e)  for  in-tree  partial 
orders  is  introduced  in  [JN83].  Here  P*  is  the  optimal  value,  and  the  log(P*)  factors  can  be  replaced  by 
log(n)  in  the  time  bounds. 

The  following  theorem  generalizes  the  results  of  [JN83]  by  allowing  multi-criteria  objective  functions. 

Theorem  12  Let  ’I'  be  the  collection  of  feasible  regions  for  the  Partially  Ordered  Knapsack  problem  with 
general  tree  partial  orders  and  let  oj  G  G”.  Then  the  problem  (w)  ,w.  Opt)  has  an  FAS. 

Proof.  The  standard  dynamic  programming  approach  can  be  used  to  solve  the  STT  Network  Flow  problem 
()$,AV  j  Opt)  in  VPP  time.  An  explanation  of  the  details  can  be  found,  for  example,  in  [JN83]. 

Applying  Lemma  2  and  Lemma  4  in  [S094]  with  Theorem  2  proves  this  theorem.  ■ 

Reversing  the  constraints  that  define  the  feasible  region  produces  a  Partially  Ordered  Knapsack  Covering 
problem  with  precedence  constraints  the  reverse  of  those  in  the  original  problem.  Application  of  Theorem  10 
would  yield  a  solution  to  the  wrong  problem. 

Note,  however,  that  if  the  original  precedence  constraints  form  a  tree  partial  order,  then  the  reversed  con¬ 
straints  also  form  a  tree  partial  order.  If  the  original  precedence  constraints  are  reversed  before  the  covering 
transformation  is  applied,  then  a  covering  problem  with  the  original  precedence  constraints  results.  The 
following  conclusion  is  therefore  valid,  even  though  not  all  the  original  constraints  are  reversed  to  pose  the 
covering  form  of  a  Partially  Ordered  Knapsack  problem. 

Corollary  4  The  Partially  Ordered  Knapsack  Covering  problem  with  general  tree  partial  orders  has  an  FAS. 
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9  Production  Planning  Problems 


The  production  planning  problem  considered  here  is  a  capacitated  single  commodity  dynamic  lot  sizing 
problem.  The  deterministic,  dynamic  demand  must  be  satisfied,  but  may  be  backlogged.  Time-dependent 
production  capacities  are  specified.  The  object  is  to  determine  how  many  units,  if  any,  of  the  item  to  produce 
each  period  in  order  to  minimize  the  (possibly  non-linear)  costs  of  production  and  inventory  storage  over  a 
finite  horizon. 


Given  demands  dj,  production  capacities  Cj,  and  production  quantities  Xj,  define  the  following  cumulative 
quantities: 

j 

Di  =  E  di  cumulative  demand  through  period  j, 

i  =  l 
3 

Ci  =  '£<“  cumulative  capacity  through  period  j, 

i-l 

j 

X,  =  cumulative  production  through  period  j,  and 


«  =  1 

Ij  =  Xj  —  Dj  inventory  at  the  end  of  period  j. 


The  feasible  set  for  this  problem  is  defined  as  follows: 


Capacitated  Single  Commodity  Dynamic  Lot  Sizing  Feasible  Set 

Instance:  (n,  d,  c,  s,  t).  The  number  of  periods  n;  the  demands  dj,  j  =  1, ..  .,n]  the  production  capacities 

Cj,  j  =  l,...,n;  the  inventory  capacities  Sj ,  j  —  l,...,n;  and  the  backordering  capacities  tj, 
j  =  1, . . . ,  n;  are  non-negative  integers  such  that  Dj  <  Cj,  j  =  1, . . .  ,n. 

Feasible  Set:  S{n,  d,  c)  is  the  set  of  all  x  G  Z”"*"  satisfying 

Xj  ^  Cj ,  j  —  1 , . . . , 

J  j  =  1, . . .  ,n 
Ij  >  ~tj ,  j  —  1, . . .  ,n 
Iq  =  0 
In  =  0 

The  assumption  that  Dj  <  Cj  is  sufficient  for  the  feasible  set  to  be  non-empty[FK71,  FLRK80].  The 
constraint  that  /q  =  0  can  be  removed,  if  desired;  simply  add  another  period  j  =  0  at  the  beginning  in  which 
xo  —  lo  is  produced  at  a  suitably  small  cost[FLRK80].  The  constraint  that  /„  =  0  is  added  without  loss  of 
generality [Zan66],  and  is  commonly  used  to  simplify  calculations. 
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The  cost  functions  typically  used  with  this  problem  are: 

Pj{xj)  cost  of  producing  xj  units  in  period  j, 

hj{Ij)  cost  of  holding  Ij  units  of  inventory  at  the  end  of  period  j,  and 
bj{Ij)  cost  of  backordering  Ij  units  in  period  j. 

This  is  an  important  problem  that  has  received  a  lot  of  attention  in  the  literature.  A  good  review  of  previous 
results  for  this  and  related  problems  can  be  found  in  [BRG87].  A  snapshot  of  the  state-of-the-art  with  regard 
to  known  complexity  results  for  variants  of  this  problem  appears  in  [FLRK80].  That  paper  has  provided  a 
framework  for  classifying  many  subsequent  results. 

The  uncapacitated  case,  i.e.  Cj  —  oo,  with  concave  functions  pj  and  hj,  was  first  treated  in  [WW58].  An 
O  (n®)  algorithm  for  the  problem  was  presented  in  that  paper.  Backlogging  for  a  limited  number  of  periods 
is  accommodated  by  an  O  (n^)  algorithm  that  appears  in  [Zan66].  An  O  (n‘^)  algorithm  for  the  capacitated 
problem  with  constant  capacities  is  described  in  [FK71];  time-dependent  capacities  are  allowed  in  [LL82]. 

Linear  production  and  holding  costs  and  time-dependent  set-up  costs  are  handled  by  an  algorithm  in 
[BDMS78];  it  runs  in  time  O  (n^D„  +  nDnCn)-  A  modification  that  allows  time- depen  dent  capacities  and 
backlogging  and  runs  in  time  O  {n^Dn  +  n^DnCn)  is  described  in  [FLRK80].  More  general  production  and 
holding  costs  are  treated  in  [BY82].  The  performance  of  various  heuristics  is  studied  in  [BMY84,  BM86], 
and  an  FAS  for  a  continuous- time  problem  appears  in  [GJ90]. 

Note  that  quantities  other  than  the  production  amount  can  be  restricted.  In  [Lov73] ,  for  instance,  a  maximum 
and  minimum  inventory  level  for  each  period  is  specified.  Also,  most  authors  assume  that  the  production 
cost  functions  are  concave,  most  often  consisting  of  a  component  consisting  of  a  set-up  cost  and  a  component 
that  is  linear  in  the  amount  produced.  An  exception  is  [Vei64],  in  which  convex  cost  functions  are  assumed. 

The  time  bounds  for  previous  algorithms  for  the  problem  considered  here,  although  pseudo-polynomial,  are 
not  VPP.  This  means  that  none  of  them  can  be  directly  converted  into  an  FAS  in  the  manner  described  in 
[S094].  The  only  exception  is  the  FAS  presented  in  [D081]. 

Nevertheless,  as  the  following  theorem  shows,  the  Capacitated  Single  Commodity  Dynamic  Lot  Sizing  prob¬ 
lem  does  have  an  FAS.  It  therefore  also  has  a  VPP  algorithm  which,  for  some  situations,  may  be  preferable 
to  the  pseudo-polynomial  algorithms  already  known. 

Theorem  13  Let  ^  be  the  collection  of  Capacitated  Single  Commodity  Dynamic  Lot  Sizing  feasible  sets  and 
w  £  fi.  Then  (w)  ,w.  Opt)  has  an  FAS. 

Proof.  Let  Hi  =  (^,  MA4  (w)  ,  w,  Feas)  and  112  =  (C2,  MA4  (w)  ,  w,  Feas),  where  the  feasible  set  C2  = 
{Al/CG/Ml/N D/FE).  The  reduction  Hi  ocvpp  112  will  be  shown.  By  Lemma  4  in  [S094]and  Theorem  7, 
this  will  prove  the  present  theorem. 

Suppose  that  an  instance  7i  =  (5, /,w,M)  of  Hi  is  given,  with  S  =  S{n,d,c).  The  objective  function  /  is 
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the  sum  over  the  decision  time  periods  of  the  functions  pj ,  hj ,  and  bj  described  above.  Define  an  instance 
I2  =  (C2,  /,  w,  M)  of  112,  where  C2  is  defined  by  an  STT  network  G  =  (A^,  A,  u,  p,  b).  The  set  N  has  (n  -f  1) 
nodes:  a  node  jSj  corresponding  to  each  period  j,  with  demand  dj;  and  node  0.  The  set  A  has  (3n  —  2)  arcs 
that  are  described  in  the  following  table: 


Arc 

Capacity 

Multiplier 

Value 

(0, /?;■)> 

i  =  1,  • 

..,n 

1 

i  =  1,  • 

. . ,  n  —  1 

Sj 

1 

i  =  1, . . 

. ,  n  —  1 

^3 

1 

_ i 

The  production  in  period  j  corresponds  to  the  value  xoj,  and  the  cost  of  a  production  schedule  is  the  same 
as  the  cost  of  the  corresponding  flow  in  C2  ■  The  reduction  associates  a  solution  of  Ii  with  a  solution  of  I2 
that  has  the  same  objective  function  value.  ■ 


10  Summary 

This  work  introduced  the  STT  Network  Flow  problem,  which  generalizes  a  wide  variety  of  commonly-studied 
combinatorial  optimization  problems.  A  VPP  algorithm  for  solving  a  simple  version  of  this  problem  was 
demonstrated  and  then  used  to  show  that  the  problem  has  an  FAS.  The  boundary  between  easy  and  hard 
problems  on  STT  networks  was  explored  in  some  detail. 

The  computational  complexity  of  finding  approximate  solutions  to  several  problems  that  arise  frequently 
in  applications  was  also  studied.  The  previously  known  FAS  for  the  binary  case  of  the  single- criterion 
Arborescent  Knapsack  problem  was  generalized  to  the  integer  case  with  multiple  criteria.  An  FAS  for  the 
covering  form  of  this  problem  was  also  demonstrated;  this  problem  has  not  been  considered  before  except 
for  the  special  case  of  the  Min  Job  Sequencing  with  Deadlines  problem. 

An  FAS  for  the  multiple  criteria  form  of  the  Partially  Ordered  Knapsack  problem  with  general  tree  constraints 
was  also  shown,  thereby  generalizing  the  previously  known  single-criterion  results.  The  covering  form  of  this 
problem,  which  has  not  been  discussed  previously,  was  also  shown  to  have  an  FAS. 

Finally,  the  Capacitated  Single  Commodity  Dynamic  Lot  Sizing  problem  was  considered.  The  previously 
known  pseudo-polynomial  time  algorithms  for  this  problem  were  not  VPP,  and  so  could  not  be  used  directly 
to  show  the  existence  of  an  FAS.  By  reducing  this  problem  to  an  STT  Network  Flow  problem,  however,  the 
existence  of  an  FAS  was  demonstrated.  This  also  implies  the  existence  of  a  VPP  algorithm  that  might  be 
more  useful  in  some  situations  than  are  previously-known  algorithms. 
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